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Abstract 

In my diploma thesis I calculate a lower mass bound of the lightest neutralino assuming 
an additional Supernova cooling mechanism. I consider elektron-positron annihilation and 
nucleon-nucleon-bremsstrahlung as cooling processes. These cooling processes must not 
be too effective, because tliey would destroy the neutrino signal (SN 1987a). From this 
requirement I derive an upper bound for the emitted energy and as a consequence I get 
a lower bound for the mass of the lightest neutralino. 
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Uberblick 



In der vorliegenden Arbeit habe ich eine untere Massenschranke für das leichteste Ncutralino 
berechnet. Bei einer Supernova- Explosion wird der größte der Teil der freiwerdenden 
Energie in Form von Neutrinos abgestrahlt. Dieses Signal ist durch Simulationen bestätigt 
worden. Bei der Supernova SN 1987 A sind diese Neutrinos zum ersten Mal nachgewiesen 
worden. Ihr Signal dauerte ungefähr 10 s. Ein zusätzlicher Kühlmechanismus darf nicht 
zu effektiv sein. Er würde die Energie, die die Neutrinos wegtragen müssen, schon vorher 
abführen. Dann würde man ein kürzeres Neutrinosignal messen. Folglich läßt sich aus der 
Forderung, dieses Signal nicht zu zerstören, eine obere Schranke für die maximale Stärke 
eines zusätzlichen Kühlmechanismus ableiten, nämlich Kühlung durch die Emission von 
Neutralinos. Dies geschieht ohne eine Supernova zu simulieren. 

Ausgangspunkt sind zwei Produktionsmechanismen für Neutralinos im Supernova- 
Core^, nämlich Vernichtung von Elektronen und Positronen sowie Nukleon-Nukleon— 
Bremsstrahlung. Das Gore wird als durchsichtig für die Neutralinos angenommen. Diese 
Annahme ist für einen Selektronenmasse im Bereich von mindestens 1 TeV erfüllt. Für 
den ersten Mechanismus ist ausgehend vom totalen Wirkungsquerschnitt die abgestrahlte 
Energie in Abhängigkeit der Masse des leichtesten Neutrahnos berechnet worden. Für die 
Temperatur und die Entartung im Gore der Supernova sind Profile aus einer gegebenen 
Simulation gewählt worden. Diese Simulation hat die Neutralinos nicht berücksichtigt. Die 
Masse des Selektrons ist zwischen 100 GeV und 1000 GeV liegend angenommen worden. 
Die Abstrahldauer der Neutralinos ist im Bereich von 0.002 — 2s variiert worden. Um 
die Massenschranke zu erhalten, wurde die abzustrahlende Energie nach oben begrenzt. 
Die Mischungsparameter des leichtesten Neutralinos haben bei kleinen Selektronenmasscn 
keinen Einfluß auf die abgestrahlte Energie. Das leichtetste Neutralino besteht zu mindestens 
98% aus Bino, der Rest verteilt sich hauptsächlich auf die beiden Higgsinos. Mit diesem 
Verfahren erhält man untere Schranken im Bereich von 50MeV (Selektronmasse groß, 
Abstrahldauer klein) bis 350 MeV (Selektronmasse klein, Abstrahldauer groß). 

Im zweiten Mechanismus ist die Emissivität direkt mit dem Raffelt-Kriterium - einer 



^Ich verwende "Gore" für das aktive Sternzentrum und "Kern" als Kurzform von Atomkern. 
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TABELLENVERZEICHNIS 



oberen Schranke für die Emissivität, die durch Simulationen überprüft worden ist - verghchen 
worden. Da die Strahlung als "sofl' genähert wird, erhält man einen Ausdruck für die 
Emissivität in Abhängigkeit der Neutralinomasse und der Coretemperatur, der für alle 
Teilchen verwendbar ist, die an den axialen Strom koppeln. Hier ist die Squarkmasse 
wie oben beim Selektron eine unbekannte Größe, für die ein ähnlicher Massenbereich wie 
für die Selektronenmasse zugrunde gelegt wird. Damit erhalte ich Massenschranken in 
derselben Größenordnung. 

Ein leichtes Neutralino im Bereich von 34MeV kann somit nicht durch Supernova- 
Kühlung ausgeschlossen werden, sofern die Sfermionmassen hoch genug sind. 



Kapitel 1 
Einführung 



1.1 Einführung in die Themenstellung 

Das Ziel meiner Arbeit ist die Berechnung einer unteren Massenschranke für das als 
leichtestes supersymmetrisches Teilchen angenommene Neutralino Xi mittels Supernova- 
Kühlung. Dabei darf die Supernova nicht zu viel Energie in Form von Neutralinos abstrahlen, 
damit das beobachtete zeitverzögerte Ncutrinosignal nicht zerstört wird. Als 1987 der 
Stern Sanduleak 69202 in der Großen Magellanschen Wolke explodierte, konnten zum 
ersten Mal die bei einer Supernova entstehenden Neutrinos in den Wasser- Tscherenkow- 
Detektoren Irvine-Michigan-Brookhaven (IMB) und Kamiokande II beobachtet werden. 
Diese waren eigentlich zur Messung des Protonzerfalls gebaut worden. Neutrinos tragen 
wesentlich zur Kühlung einer Supernova-Explosion bei. Sie führen hundert bis tausend 
Mal so viel Energie wie Photonen ab. Ihr Signal dauerte etwa 10 sec. Jeder zusätzliche 
Kühlmechanismus darf nicht zu effektiv sein. Denn sonst würde dieser Kühlmechanismus 
das Gore so weit kühlen, daß die Neutrinos nicht mehr gefangen wären und das Neutrinosignal 
kürzer wäre. Das steht jedoch im Widerspruch zu den Beobachtungen. Im folgenden 
nehmen wir an, daß es einen weiteren Kühlmechanismus gibt, der Energie aus dem 
Supernova- Gore nach außen abführen kann, nämlich die Produktion von Neutralinos. 
Neutralinos können auf folgende Arten erzeugt werden: 

• Elektron-Positron-Annihilation in ein Z'^-Boson: e~^e~ — > — XiXi 
Elektron-Positron-Streuung: e+e" XiXi unter Austausch eines Selektrons 

• Nukleon-Nukleon-Bremsstrahlung: NN NNxiXv 

Bei diesem Prozeß wird entweder ein Z*^ abgestrahlt, das in zwei Neutralinos zerfällt, 
oder ein Quark geht kurzzeitig in ein Squark über und strahlt dabei zwei Neutralinos 
ab. 
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Für diese Mechanismen berechnen wir die Emissivität (Energieabgabe pro Zeit und Volumen 
oder Energieabgabe pro Zeit und Masse) in Abhängigkeit von der Neutrahnomasse. Im 
ersten Fall liefern Volumen- und Zeit Integration die absolute freigesetzte Energie in Abhängigkeit 
von der Neutrahnomasse. Diese Berechnungen erfolgen numerisch. Man kann aus der 
Anfangs- und der Endgröße des Gores abschätzen, wieviel Energie von einer Supernova 
maximal abgestrahlt wird. Dies sind etwa 3-10^^ erg. Den Großteil davon führen Neutrinos 
ab. Im Durchschnitt führt jede der sechs Neutrinoarten demzufolge 5 ■ 10^^ erg ab. Die 
Annahme, daß 10^^ erg in Form von Neutralinos abgestrahlt werden, sollte nicht zu 
einer Zerstörung des Neutrinosignals führen. In die Berechnungen gehen eine Reihe von 
Parametern ein: Temperatur, chemisches Potential und Dichte des Sterncores, die Mischungsparamete: 
des leichtesten Neutralinos, die Selektronmasse und die Dauer der Energieabstrahlung. 
Für diese Parameter wurden geeignete Modelle oder Werte zugrunde gelegt. 

Im Falle der Nukleon-Bremsstrahlung ist es nicht notwendig, die Emissivität e über das 
Volumen und die Zeit zu integrieren. Hier kann die Emissivität sofort mit einem Kriterium 
von Raffelt verglichen werden: e < 10^^ erg/ (g-s) (aus |Raf96j ) . Dieses Kriterium ist durch 



Simulationen bestätigt worden. Das Endergebnis ist für andere Theorien verwendbar. 



1.2 Was sind Neutralinos? 
1.2.1 Anschauliche Beschreibung 

Zuerst soll kurz erläutert werden, was Neutralinos überhaupt sind. Supersymmterie ist 

eine Symmetrie zwischen Fermionen und Bosonen. Damit kann das Hierarchieproblem 

in der Teilchenphysik gelöst werden(s. |Mar97j . |Hab93j ) . In der Supersymmetrie (kurz 

SUSY) wird jedem Boson ein Fermion als Superpartner zugeordnet und umgekehrt. Dadurch 

wird die Anzahl der Teilchen verdoppelt. Die Eichbosonen der Elektroschwachen Kraft 

sind nach spontaner Symmetriebrechung durch ein komplexes Higgs-S'f/(2)-Dubletts die 

geladenen VT-Bosonen, das neutrale Z°-Boson und das ebenfalls neutrale Photon (s. 

|HM84p . Von den vier Freiheitsgraden dieses Higgsdubletts werden drei durch den Higgsmechanismus 

"aufgegessen" . Ein Higgs- Teilchen bleibt übrig. Zu den vier Austauschteilchen gibt es im 

"Minimalen Supersymmetrischen Standardmodell" (MSSM) vier fermionische Gauginos 

als Superpartner: die geladenen Winos und die neutralen Teilchen Zino und Photino (s. 

|HK85j ). In der Supersymmetrie muß man zwei komplexe Higgs-S'f/(2)-Dubletts einführen, 

damit u. a. sowohl up- als auch down-Quarks Masse erhalten f [HK85j . [GH86j ). Es bleiben 

am Ende nur fünf reelle Felder übrig, da drei Felder "aufgegessen" werden, um den 

schweren Vektorbosonen Masse zu geben. Von diesen fünf Feldern sind zwei geladen 

und drei neutral. Die fermionischen Superpartner dieser Higgs-Felder heißen (nacd der 

Brechung von SU{2)l x f/(l)y — ^ U{l)em) Higgsinos. Die zwei GP-geraden neutralen 
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Higgsinos mischen mit den neutralen Gauginos Photino und Zino sowie die geladenen 
Higgsinos mit den Winos. Die Eigenzustände der Schwachen Wechselwirkung sind also 
nicht die propagierenden Masseneigenzustände. Die neutralen Masseneigenzustände werden 
Neutralinos genannt, die geladenen Charginos. Um die Masseneigenzustände zu erhalten, 
muß man die Massenmatrix diagonalisieren. 

Des weiteren gehe ich von erhaltener i?-Parität aus (weswegen die Neutrinos nicht mit 
den Neutralinos mischen) und nehme an, daß bei hohen Energien keine GUT-Symmetrie 
auftritt, also die SU{3)c x SU (2) x f/(l)y-Theorie nicht zur 5'f/(5)-Theorie vereinigt wird. 

1.2.2 Mathematische Beschreibung 
Neutralinos 

Für die Einführung in die SUSY- Theorie verweise ich auf die Literatur: Der Superfeld- 
Formalismus, mit dem man supersymmetrische Lagrangdichten konstruieren kann, wird 
z. B. in |BL94j dargestellt. Die Konstruktion supersymmetrischer Lagrangefunktionen für 
eine SU{2) l x ?7(l)y-Theorie wird in ^HK85j beschrieben, der Higgs-Sektor wird in |GH86j 
erläutert. 

Ich fasse hier die Resultate über die Neutralino-Massenmatrix kurz zusammen: 



• Massenterme: In der Lagrangedichte des MSSM findet man folgende Massenterme 
für die Superpartner der neutralen Eichbosonen und neutralen Higgsfelder unter 
Berücksichtigung des Symmetriebruchs von SU {2) x U{1)y zu [/(l)em und der soß 
susy-breaking Terme (siehe dazu jHK85j ) (in Weylspinor-Darstellung): 

J^M = h^/{^^T^\zM]j^-V2ijjiJ + ^{Mcos^9w + M'sm^9w)\zXz 

+ (M - M') sin Ow cos ^v^AzA^ + ^(M' cos^ Ow + M sin^ 0,^) A^A^, + fi^n^^l^ 
+ h.c. . (1.1) 

Hier bedeuten g, g': Kopplungskonstanten der SU{2) x [/(l)y-Theorie, Wi, V2 sind 
die Vakuumserwartungswerte der zwei neutralen CP-geraden Higgsfelder, 6w ist der 
elektroschwache Mischungswinkel, M, M', /x sind Wino- und Binomasse sowie der 
Higgsmischungsparameter aus dem Superpotential und A^, A^, V'j/i! ^1^2 bezeichnen 
Photino-, Zino- und die beiden Higgsinof eider. 

• Basis: Anstatt der Photinobasis {A^, A^, '?/')/^, V'l^j} "wird häufig auch die Bino-Basis 
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{X' , X^,il)jj_^,ip'jj^} verwendet. Dann sieht Gleichung so aus: 

+ ^MX^X' + ^MX'X' + ^iiljh^K^+^-C- ■ (1-2) 

Mischungsmatrix: Schreibt man die Elemente der Bino-Basis als tpj, j = 1 . . .4, so 
kann man die letzte Gleichung in die Form 

= lii^^)Y.j^^ (1.3) 

bringen (über doppelt auftretende Indizes ist zu summieren) mit der (im allgemeinen 
komplexwertigen) Mischungsmatrix Y: 

( M' —mz cos ß sin 6w niz sin ß sin 6w \ 

M mz cos ß cos 6w —mz sin ß cos 6w 

-mz cos ß sin 6w mz cos ß cos 6w —ji 

\ mz sin ß sin 6w —mz sin ß cos Oy/ — /i J 

(1.4) 

tan /3 = ^ ist das Verhältnis der Vakuumserwartungswerte fi,f2 der Higgsfelder. 
M, M', /i, tan /5 sind unabhängige Parameter im MSSM. 

Diagonalisierung der Mischungsmatrix: Die Mischungsmatrix ()1.4j) ist nicht diagonal. 
Um die propagierenden Masseneigenzustände x^i erhalten, muß man die Mischungsmatrix 
diagonalisieren. Dabei können die Eigenwerte negativ werden. Deshalb nimmt man 
eine unitäre Diagonalisierungsmatrix A^, wobei man die freien Phasen so wählt, daß 
die Eigenwerte der Diagonalmatrix Nj:, positiv sind (siehe |MP90j ). Dann ist eine 
Interpretation der Eigenwerte als Massen möglich ist ([HK85]). Die Vorzeichen der 
Massen können als CP-Eigenwerte gedeutet werden f J,BNP84j ). Es gilt: 

Yd = N*YN-\ (1.5) 
X° = Nijip^, i,j = 1...4: Neutralino-Masseneigenzustände . (1.6) 

Der Massenterm ll . 3l läßt sich dann als — ^ J2i MjXi^x^schreiben, wobei Mi die Eigenwerte 
der Matrix ()1.4|) sind. Man kann dies auch in vierkomponentige Majorana-Spinoren 
umschreiben: 

Die Eigenwerte der Massenmatrix bestimmt man am besten numerisch. Es gibt aber 
auch analytische Lösungen, siehe dazu |KAS92j . 
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Charginos 

Die Mischungsmatrix der Charginos hängt von M, /i und ß ab. Weil die Gauginos 
elektromagnetisch wechselwirken können und sie deshalb im Prinzip einfacher zu entdecken 
sind als die elektrisch neutralen Neutralinos, müssen die Parameter der Mischungsmatrix 
so gewählt werden, daß die Charginos schwerer als 120 GeV sind, der momentanen unteren 
Massenschranke, die durch die Beschleuniger gesetzt ist. Hier eine Zusammenfassung der 
Resultate für die Charginos: 

• Der Chargino-Massenterm (Weyl-Spinor-Schreibweise) ist gegeben durch: 

Cm = ^ {viX+^jj, + v^X-^h,) + M\+\- - + h.c. (1.8) 

• Gl. dLH) läßt sich durch eine Basistransformation umformen: 

= (-iA+.^iJ, <;'- = (-iA-,'(/4,) (1.9) 

= ■r)(::^-)+h.c., (1.10) 

wobei die Massenmatrix X gegeben ist durch: 

X = ( ^ m^V2sin/?Y ^ ^l) 

\m^-\/2 cos ß n J ^ 

ml = ]g\vl + vl); tan/5 =^ (1.12) 
4 Vi 



Die Masseneigenzustände erhält man durch Diagonalisieren der Matrix 



X^ 
X 



Xt = V,,^l, x^ = U,,^T, CM = ~{Xi{XDhxt + ^-C-) (1-13) 
mit U, V so, daß U*XV^^ = Xc diagonal 

Aus den Weylspinoren lassen sich Diracspinoren xi, X2 konstruieren: 

XI = 1 > X2= 1 , Cm = -(M+^iXi + M_J,X2) ■ (1.14) 
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Anmerkung: Die oberen Indizes +, — an den Spinoren kennzeichnen in diesem 
Abschnitt die Ladung und nicht die hermitesche Konjugation. 

Anmerkung zu Majorana-Massen: 

Massenterme für Diracspinoren haben die Form mipip — mlipj^ipR+h. c), (■0: Diracspinor, 
'4^L,R = Pl,r'^, Pl,r = I (1 =F 75))- V'i uiid V'-R sind unabhängig voneinander und transformieren 
sich unter Lorentztransformationen verschieden. Solche Dirac-Massenterme mischen hnks- 
und rechtshändige Felder. 

— T 

Majoranaspinoren haben die definierende Eigenschaft ^E''^ = = ^, C: Ladungskonjugationsc 
Damit gilt — und Massenterme haben die Form ^(V'lV'l + h. c). Es werden 
Komponenten gleicher Händigkeit gemischt. Dadurch ergibt sich ein anderes Transformationverhalten 
für Majoranaspinoren, da allgemeine Massenterme die Gestalt i/j'^Mip + h. c. haben, wobei 
ijj = {ipi . . . ipn) Weyl-Spinoren sind. Transformiert man Xi = ^ij'^jj wobei unitäre 
Matrix ist, dann gilt {N~^x)'^M^~^X = ^* ^N^^Xi was das auf den ersten Blick 
ungewöhnliche Transformationverhalten der Neutralinomassenmatrix erklärt. Denn um 
eine hermitesche Matrix A mittels einer unitären Transformation U auf Diagonalgestalt 
Ad TM bringen, erwartet man das Transformationsgesetz Ad — UAU^. 



Kapitel 2 
Sternentwicklung 



2.1 Bildung einer Gaswolke 

In den folgenden Abschnitten werde ich die Entwicklung von Sternen aus einer Wasserstoffwolke 
kurz dargestellen, um das mögliche Supernova-Ereignis in der Entwicklung eines Sterns 
verstehen zu können. Dabei werden die thermodynamischen Zustandsgrößen zueinander 
in Beziehung gesetzt und erklärt, wie sie die Sternentwicklung steuern. Ich stütze mich 

auf [nsH3. 

2.1.1 Gravitation gegen kinetische Energie 

Um den Kollaps einer Wasserstoffwolke geringer Dichte zu verstehen, betrachte man die 
Energieverhältnisse in einer solchen Wolke: 

E = Ekin + i^pot = -\^ = r-"^^' ~ f^G-n'p'R' (2.1) 

(M, i?, p: Masse, Radius und Dichte der Gaswolke, v"^: mittlere quadratische Geschwindigkeit 
der Wasserstoffatome, G: Gravitationskonstante). 

In einem gebundenem System muß die Gesamtenergie E negativ sein, das heißt, 
die potentielle Energie muß betragsmäßig größer als die kinetische Energie sein. Da die 
potentielle Energie stärker mit dem Radius als die kinetische wächst, kann die Wolke einen 
gebundenen Zustand erreichen. Das kann auf zwei Wegen geschehen: durch Kompression 
und durch Akkretion. Kompression bedeutet die Erzeugung von Dichtewellen, die durch 
die Wolke laufen. Diese werden zum Beispiel durch Supernova-Explosionen oder durch 
vorbeiziehende Sterne erzeugt. 
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Akkretion heißt, daß die Wolke Materie aus ihrer Umgebung einsammelt und so sich 
ihr Durchmesser vergrößert. Dies geschieht in kosmologischen Zeiträumen, das heißt, der 
Vorgang benötigt etwa eine Milliarde Jahre. Den Radius dieser kollabierenden Wolke kann 
man leicht abschätzen, indem man E zu Null setzt. Dann folgt für den kritischen Radius 

D _ 6 GM 

Der Kollaps wird durch Gegenkräfte angehalten, die im folgenden diskutiert werden. 



2.1.2 Quasihydrodynamisches Gleichgewicht 

Wir nehmen eine sphärische Massenverteilung an und betrachten eine Schale der infinitesimalen^ 
Dicke Ar vom Radius r. Deren Volumen ist Anr^Ar, ihre Masse AAfs = Anr'^pAr. Auf 
diese Schale wirkt die Kraft AFs = —QMa^Ms _ Gesamtkraft lautet demnach: 

^2 

AFges = Anr^AP = OgesAM, . (2.2) 

Daraus ergibt sich die Gesamtbeschleunigung Oges zu 



1 AP GM , ^ „ „ „ 
ages = - 7^ ^obei AP^P,-P, . (2.3) 



Hier können wir jetzt einige Fälle betrachten: 



1. ages = 0: Gleichgewicht 

2. Odes — 0: Kleine Beschleunigung, sog. quasihydrostatisches Gleichgewicht QHE 
(quasihydrostatic equilibrium) 

3. Oges > 0: expandierende Wolke 

4. Oges < 0: 

(a) ^ > 0: nicht physikalisch 

< l^^l: schlimmstenfalls freier Fall, Kollaps 



(b) 



1 AP 
p Ar 



Ein Stern (Materieansammlung) ist immer bestrebt, ins QHE zu gelangen, insbesondere 
ist eine Kontraktionsphase dann beendet, wenn der Stern sich wieder im QHE befindet. 



Aa; bezeichnet die infinitesimale Größe zu x 



2.1. BILDUNG EINER GASWOLKE 
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2.1.3 Die Gasgleichung 

Um den Tatsachen Rechnung zu tragen, daß ein Stern erstens nicht nur aus einer Teilchensorte 
besteht, zweitens die Teilchen manchmal ionisiert sind und daß drittens das Volumen keine 
dem Problem angemessene Größe ist, wird die Gasgleichung PV = NkT in eine passende 
Form gebracht. Es seien x, y, z die Massenanteile H, He und Metall (= alles, was schwerer 
als He ist), also x — M-n/M, Mh: Gesamtmasse des Wasserstoffs usw. Dann kann man 
die Anzahl der Teilchen ausrechnen (nicht ionisiert): 



A^H = — = (2.4) 

Mne MHe yM 

Nue = = = (2.6) 

ruMe ArriH ruAn 

/ 1 1 \ M , , 

\ 4 A / mn 

1 



Im Falle vollständiger Ionisierung liefert ein Atom mit A Nukleonen Z Elektronen, damit 
haben wir Z+1 Teilchen. Daher müssen die vorhergehenden Gleichungen um entsprechende 
Vorfaktoren verändert werden: 

M 

iVn = 2x , (2.8) 

Nne = ly—, (2.9) 
4 mn 

Z+1 MIM , , 

A^Me = —^z 2.10 

A mn 2 mn 

3 1 , M / 1 \ M , , 

Ntot = {2x+-y+-z) = - . (2.11) 



4 2 mn V^/inn "^h 
Jetzt kann die Gasgleichung umgeschrieben werden: 



PV=l^kT. mit ,= ^^^ P = f^ (2.12) 

w mn V wmn 
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Aus dieser Gleichung kann man schließen, daß ein Druckgradient drei Ursachen haben 
kann: 

Temperaturgradient (2-13) 

Dichtegradient (2.14) 

Teilchengradient (2-15) 

Später werden wir sehen, wie all diese Gradienten dazu beitragen, daß ein Stern nicht 
kollabiert. 

2.2 Thermonukleare Prozesse 

2.2.1 Eine weitere Energiequelle: Kernfusion 

Sterne wie die Sonne beziehen ihre Energie aus der Verschmelzung von Wasserstoff zu 
Helium oder ähnlichen Prozessen. Aus dem Verlauf der Kurve der mittleren Bindungsenergie 
pro Nukleon in Abhängigkeit von der Massenzahl kann man folgern, daß man durch Fusion 
von Nukliden, die leichter als die eisenartigen Metalle sind, Energie gewinnen kann. Für 
den Nettoprozeß 

Qp — >'^Re + 2p + 2e+ +Ü+ E (2.16) 

sind Temperaturen von 10^ K ^nötig. Gravitation und Kernfusion sind jetzt zwei gegeneinander 
wirkende Mechanismen: Durch die Fusion erhöht sich die Sterntemperatur im Gore und 
infolge dessen auch der Druck. Dieser Druck wirkt der Gravitation entgegen. 

2.2.2 Entwicklung zum Roten Riesen 
Prozesse im Inneren eines Sterns 

Der Stern besteht aus zwei Teilen: einem Gore, in dem Wasserstoff zu Helium fusioniert 
wird, und einer äußeren, "kälteren" Hülle. Im Gore reichert sich mit der Zeit immer mehr 
Helium an, dabei wird die Teilchenanzahl kleiner, nach der Gasgleichung nimmt damit 
auch der Druck ab, und der Kern zieht sich zusammen. 

^Die Geschwindigkeit der Protonen unterliegt einer Maxwellverteilung mit einer mittleren 
Geschwindigkeit von lO'^ K. Eine solche energie würde allein nicht ausreichen, um die Coulomb- Abstoßung 
zu überwinden, aber die Anzahl der Teilchen im Schwanz der Verteilung mit einer entsprechenden 
Geschwindigkeit ist groß genug, um die Fusion aufrecht zu erhalten. 



AP AT 

Ar Ar 

AP Ap 

-IT- ^ -TT- 

Ar Ar 

AP Avj-^ 

-r — (X — 

Ar Ar 
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An der Grenzfläche stimmen Druck und Temperatur von Hülle und Gore überein, aus 
der Gasgleichung folgt dann ^ = Da Wi mit der Zeit wächst, nimmt die Dichte des 
Gore zu. Seine Masse bleibt aber konstant (den Massendefekt kann man vernachlässigen), 
so daß das Volumen des Sterns schrumpft. Bei der Kontraktion des Gores nimmt jedoch die 
Gravitationsenergie weiter ab, die freigesetzte Gravitationsenergie erhöht die Gesamtleuchtkraft 
des Sterns, und der Temperaturgradient im Reststern wird größer. Aus Gl. ()2.3|1 und der 
Gasgleichung folgt, daß die Hülle sich ausdehnen muß: Wir haben ein kontrahierendes 
Gore und eine expandierende Hülle. Der Stern kann sich dabei nicht aufheizen, da die 
Energie zur Expansion verwendet wird. 

Es gibt aber noch einen anderen Prozeß. An der Oberfläche des Sterns wächst mit 
der Expansion die potentielle Energie eines Teilchen. Weil die Energie aber konstant ist, 
muß seine kinetische Energie abnehmen, infolgedessen auch seine Oberflächentemperatur. 
Währenddessen nimmt im Gore die Anzahl der Protonen ab. Da sich die Emissivit ät wie 
Spp = x^T^ verhält, nimmt die Fusionsrate ab {x wird kleiner) damit einher geht eine 
verminderte Kerntemperatur und Leuchtkraft. Dies wiederum führt zu einer Kontraktion 
des Sterns, wodurch sich ein neues Gleichgewicht einstellen kann. Es gibt also zwei gegeneinander 
wirkende Prozesse, welcher davon die Oberhand gewinnen kann, hängt von der Masse des 
Sterns ab. Die weitere Entwicklung des Sterns hängt von seiner Masse ab. 



Alterung des Sterns 

Die Alterung beginnt mit der Bildung eines festen homogenen Heliumcores. Die Temperatur 
ist für ein Zünden des nächsten möglichen Prozesses, nämlich 3a ^^G, zu niedrig (10^ 
K^ sind nötig, a = ^He^^). Das Gore ist somit tot, seine Leuchtkraft Null, man sagt, 
daß das Gore jetzt isotherm sei (da ^ = 0). Der Stern würde gerne im QHE bleiben, 
im Gore gibt es keinen Temperaturgradienten, also auch keinen daraus resultierenden 
Druckgradienten. Aber der Stern kann einen Dichtegradienten entwickeln: 

dP dp kT GM 
dr dr wmu 

Dies geschieht so: Das isotherme Gore kann die Hülle nicht länger tragen, der Stern 
kontrahiert, dabei wird Gravitationsenergie frei, und es besteht die Möglichkeit, daß 
dadurch eine Schale um das Gore herum eine Temperatur zum Fusionieren von Wasserstoff 
erreicht. 

Diese Schale sorgt jetzt für den Dichtegradienten und die Entwicklung zum Roten 
Riesen: Sollte der Radius der Schale schrumpfen, bewegt sie sich in Richtung größerer 
Dichte und Temperatur. Da die Emissivität proportional mit der Dichte wächst, wird die 

siehe eben 
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Energieproduktion so lange größer, bis Energieabgabe (oc L) und Energieproduktion (oc 
Epp) sich ausgleichen. Bei der Gleichgewichtstemperatur ist der Druck nach außen größer, 
so daß die nach innen gerichtete Kraft kompensiert werden kann. Eine Störung dieser 
Gleichgewichtslage führt zu einer gedämpften Schwingung um diese Position und deshalb 
bleibt der Schalenradius konstant. Da rs sich nicht ändert, ist auch das eingeschlossene 
Corevolumen und die Coremasse konstant, folglich auch die mittlere Dichte. Der Strahlungsstrom 
aus der Schale richtung Gore drückt diesen zusammen. Dadurch wird seine Dichte größer, 
die Dichte der Schale muß kleiner werden. Es gibt deshalb einen Dichtegradienten, der 
den Einsturz des Sterns verhindern kann. 

Aber es passiert noch mehr: Da die Masse der Hülle sich ebenfalls nicht ändert, die 
Dichte an der Oberfläche Null ist und die mittlere Dichte abnimmt, dehnt der Stern sich 
aus, solange das Gore zusammengedrückt wird. Damit steuert die Schale zugleich das Gore 
und die Hülle. Durch das Ausdehnen der Hülle wird die Oberfläche des Sterns kühler. Ab 
einem bestimmten Punkt wird die Energiezufuhr in die Hülle von Konvektion bestimmt. 
Das stabilisiert die Oberflächentemperatur. Der Stern erscheint rötlich, er ist zu einem 
Roten Riesen geworden. 

2.2.3 Entartung der Elektronen 
Entctrtung der Materie 

Entartung heißt, daß die mittlere kinetische Energie der Fcrmionen - in unserem Falle 
sind das die Elektronen - kleiner als die Fcrmi-Energie ist. Fast alle Energiezustände 
unterhalb der Permi- Energie sind besetzt. Aus der Temperatur und der Teilchendichte 
läßt sich abschätzen, ob ein Gas entartet ist oder nicht. Fermi-Impuls und Fermi-Energie 
sind gegeben durch: 



Betrachten wir als Beispiel die Sonne. Ihre mittlere Dichte beträgt p ^ 1.4- lO^kg/m^. 
Ihre Masse wird hauptsächlich von Protonen getragen, deren Masse 1.67-10"^^ kg beträgt. 
Somit befinden sich 10^° Protonen in einem m^, was auch der Elektronendichte entspricht, 
da die Sonne elektrisch neutral ist. Die Fusion findet bei 10^ K = 800 eV statt, man 




(nichtrelativistisch) 



(relativistisch) 



(2.18) 



(2.20) 



(2.19) 
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kann also nichtrelativistisch rechnen. Der Fermi-Impuls beträgt ~ 6200 eV, und damit 
die Fermi-Energie 38 eV für die Elektronen und 2 ■ 10~^ eV für die Protonen. Beide Gase 
sind also nicht entartet. Beim Ubergang eines Sterns wie der Sonne zum Roten Riesen 
kann der Fermi-Druck einen beträchtlichen Beitrag zum Gesamtdruck leisten: Wir nehmen 
einen Coreradius von jqRq = 10'^ m an, eine Masse von 0.2 und eine Temperatur von 
2 • 10^ K = 2 ■ 10^ eV. Daraus folgt eine Teilchendichte (Helium) von ?7,He = 16 ■ 10^^ m~^ 

1 /3 

und ne = 32 ■ 10^^ m~^. Der Fermi-Impuls ist pF,e = ch (^) ~ 10^ eV, woraus eine 
Fermi-Energie von 10^ eV folgt. Folglich sind die Elektronen teilweise entartet. 

Es hängt jetzt von der Masse des Sterns ab, ob der Fermi-Druck der Elektronen (oder 
der Neutronen) den Kollaps des Sterns verhindern kann. Wenn der Stern schrumpft, 
werden die Elektronen zum Atomkern hin gedrückt, d. h., der ihnen zur Verfügung 
stehende Raum wird kleiner. Nach Gl. ()2.18p und ()2.19|) wächst ihre kinetische Energie 
und damit auch der Fermi-Druck. 

2.2.4 Weitere Entwicklung 

Gemäß Arnett |Arn96j kann man sich die weitere Entwicklung des Heliumcores näherungsweise 
so wie die Entwicklung des Wasserstoff fusionierenden Sterns vorstellen. Der nächste 
Prozeß nach der Fusion von Wasserstoff ist der 3a-Prozeß: 

3a — > + 7 (2.21) 
^^C + a — > (2.22) 

Aber zuerst muß die Entartung des Gores überwunden werden. Wie oben erwähnt besteht 
das Gore eines gealterten H-fusionierenden Sterns aus festem Helium. Das Gore wird aber 
von den darüber liegenden Materieschichten zusammengedrückt und heizt sich dabei auf. 
Es kann sich aber nicht wie ein Gas ausdehnen. Sobald eine Temperatur von 1.5 ■ 10^ 
K erreicht ist, setzt der 3a-Prozeß ein. Die zusätzliche Energie kann nicht schnell genug 
abgeführt werden, obwohl die Wörmeleitfähigkeit eines entarteten Elektronengases gut 
ist: Es kommt zu einer Explosion, dem sogenannten Helium-Flash. Dadurch wird die 
Goredichte und die Temperatur vermindert, die Elektronen sind nicht mehr entartet. 
Die Gravitation ist aber so stark, daß sie die Explosion so weit dämpfen kann, daß kein 
Material ausgeschleudert wird. Es setzt eine weitere Phase der Kontraktion ein. Dadurch 
kann im Gore des Sterns, der zur Zeit keine nukleare Energiequelle besitzt, wieder die 
Temperatur für den 3a-Prozeß erreicht werden. Die Kernfusion läuft diesmal kontrolliert 
ab, weil das Gore gasförmig und nicht fest ist. Die bei der Wasserstoffusion entstandene 
Asche ist also in einem Folgeprozeß zur Energiegewinnung genutzt worden. 

Insgesamt gibt es sechs solcher Fusionsprozesse. Dies führt dazu, daß ein ausgebrannter 
Stern einen schalenförmigen Aufbau (ähnlich einer Zwiebel) besitzt, weil die Asche des 
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vorausgegangenen Prozesses der Brennstoff des nächsten ist. Das funktioniert solange, bis 
die mittlere Bindungsenergie pro Nukleon ein Maximum erreicht hat. Dieser Fall tritt bei 
den eisenartigen Metallen (Ni, Co, Fe) ein. Die Hauptprozesse sehen im einzelnen so aus 
(siehe |Arn96j . |BB85l ): 

1. 4p — > ^He + 2e+ + 2z7 bei T ^ 10^ K 

2. 3« — > + ^/2Q ^ ^ — ^ 160 bei T ^ 1, 5 ■ 10*^ K 

3. + — > 20Ne + a, oder ^^Mg^^O + ^He — > ^o^e bei T ^ 1 ■ 10^ K 

4. 20Ne + a — > ^^Mg+ bei T ^ 1, 5 ■ 10^ K 

5. + — > 32g Q^g^ 28gj + a hei T ^2-10^ K 

6. 28Si + 28si — , 56j^j — , 56pg ^ 2e+ + 2z7 bei T ^ 3 ■ 10^ K 

In Abb. 12.11 ist der schalenförmige Aufbau eines Sterns von 25 M© skizziert, der alle 
Fusionsprozesse durchlaufen hat. 

Diese Reaktionen sind nur die wesentlichen Prozesse, daneben laufen noch viele weitere 
Reaktionen ab, die für die Nukleosynthese wichtig sind. Das setzt aber voraus, daß der 
Stern ausreichend Masse hat, so daß die Gravitation den Stern genügend zusammenquetschen 
kann. Der innerste Prozeß erzeugt Eisen, das Element mit der höchsten Bindungsenergie 
pro Nukleon. Die weitere Entwicklung hängt jetzt von der Masse und der Elektronenanzahl 
pro Nukleon dieses Eisenkerns ab. Die Grenzmasse zwischen Stabilität und Instabilität 
ist gerade die Chandrasekhar-Grenze, wird sie überschritten, ist Stabilität nicht möglich. 
Wenn der Stern es schafft, in einer Supernova- Explosion genügend Masse abzustoßen, 
bleibt ein kleiner, extrem dichter Neutronenstern übrig (für diesen gilt aber ein analoges 
Stabilitätsargument wie für den Weißen Zwerg), ansonsten ein Schwarzes Loch. 



2.2.5 Die Chandrasekhar-Grenze 

Im letzten Abschnitt habe ich die Chandrasekhar-Grenze erwähnt. Nun soll an dieser 
Stelle ein Ausdruck für diese Grenzmasse hergeleitet werden. Dabei stütze ich mich auf 
|Wei72j . 



Differentialgleichung für die Sternstruktur 

Wir nehmen an, die Sternmaterie verhalte sich wie eine ideale, isotrope Flüssigkeit, d. h. 
keine Raumrichtung ist ausgezeichnet. Die mittlere freie Weglänge im Inneren des Sterns 
sei kleiner als die Längenskala eines Beobachters, d. h. die innere Struktur der Flüssigkeit 
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spiele keine Rolle. Außerdem sei die Flüssigkeit für den Beobachter in Ruhe. Dann können 
wir den Energie-Impuls- Tensor hinschreiben: 



T 



00 



u, T 



Oi 



T 



iO 



= pßij 



(2.23) 



wobei u die Energiedichte und P den Druck bezeichnen. Die zugrundeliegende Metrik soll 
statisch und isotrop sein (r: Eigenzeit): 



dr^ = F{r)dt^ - 2E{r)dtrdr - D{r)r'^dr'^ - C{r){Ar'^ + sin^ öd(^2 ^ ^2^q2 



[2.24) 



Durch eine geeignete Transformation der Zeit kann man den in der Raumzeit gemischten 
Anteil zum Verschwinden bringen: t ^ t + $(r), außerdem transformiert man r i— 
^yC{r)r. Man erhält nach einer kurzen Rechnung: 

dr^ = B{r)dt^ ~ A{r)dr^ -r\sin^edip^ + de^) mit 



A{r) 



1 



G(r) 
C{r) 



1 + 



lrC'{r) 
2 C{r) 



und der Nebenbedingung $'(r) 



B{r) = F{r) 
rE(r) 



F{r) 



(2.25) 
(2.26) 

(2.27) 



Hier ist /' die Abkürzung für Damit sieht unsere Metrik wie folgt aus: 



9ßu 



(-B{r) 
A{r) 





V 








\ 




sin^ 9 ) 



(- 



9 



B{t) 







\ 



A(r) 












\ 





1 / 



^2.28) 



sin 

Das Ziel ist es, die Einstein-Gleichungen für ein solches System zu lösen. Allgemein lauten 
sie: 



R 



1 

1^9 ßvR 



-S-kGT, 



(2.29) 



Der Energie- Impuls- Tensor auf der rechten Seite der Gleichung ist bekannt, siehe Gleichung 
()2.23|) . Ricci- Tensor und -Skalar ergeben sich aus: 



R 



R 



R^ßXu - 
9'"'Rßu 



ßX 



~ ^ ßX^ Up ^ ßu^ Xp j 



(2.30) 
(2.31) 



Die Christoffelsymbole Tx'^ ^ sind wie folgt definiert: 



^x^ß = 9"" {9x9 ßu + dß9xu - du9ßx) 



(2.32) 
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Nach einer längeren Rechnung folgen dann die Einstein-Gleichungen: 

B" B' ( N B'\ A' ^, , . , , 

R' / A' R'\ R' 

- -21 + 21 (ä + ¥ ) - M - -'^"^^ - ^^)^ '^■^^) 

-R<^(/) = -Ree (2.36) 

Rf,^ = 0, falls fi^u (2.37) 

mit der Dichte p. Zusätzlich benötigen wir eine Gleichung für hydrodynamisches Gleichgewicht. 

Aus dßT'''^ = folgt 



B' 2P' 



(2.38) 



B P + p 

Bildet man + ^ + so erhält man die Gleichung: 

r 

/r\' / 2G f \ -1 

(^-j = 1 - SnGpr^ ^A=(l-— Anx^pdxj . (2.39) 

^ ^ 

M{r) 

Dies setzen wir in Gleichung ()2.34p ein, und es folgt nach einigen Umformungen 



Newtonsche Sterne 



Wir betrachten jetzt Sterne, deren innere Energie und Druck viel kleiner als ihre Restmassendichte 
sind, wo relativistische Effekte keine Rolle spielen. Dann vereinfacht sich Gleichung ()2.4Ü|) 
zu (siehe |Wei72p : 

1 2 J 

^ P(r) = -AnGr^pir) . (2.41) 



dr p{r) dr 
Eine Variablentransformation 



p = p(0)ei/(^-i), p = Kp{QrQ^/^^-'\ r = p(0)(^-2)/^a: (2.42) 
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führt zur sogenannten Lane-Emden-Gleichung: 

4-^a;2Ae + e^/(^-^),e(o) = 1, e'(o) = o. (2.43) 

ax ax 

Für 7 > I hat die Lösungsfunktion eine Nullstelle bei xi, d. h. Druck und Dichte 
verschwinden bei xi. Damit folgen der Sternradius R und die Sternmasse M zu (mit 
den Bezeichnungen aus |Wei72p 



M = AnpiOf^-'^/'f ^ X '\-,^''-''^'^^-'\l\Q\xr)\^ . (2.45) 

V47rG(7 - 1)/ 

Bedingung für Stabilität 

Nun können wir eine Bedingung angeben, wann ein Stern stabil ist. Dazu betrachten wir 
seine Energie 

E = T + V= ^ p'-' --(-A GM''/^ _ (2.46) 

7-1 5 V3 y 



Damit ein Stern instabil werden kann, muß er in Bezug auf seine zentrale Dichte bei 
konstanter Entropie pro Nukleon und konstanter chemischer Zusammensetzung gemäß 
Theorem 1 aus |Wei72j stationär werden, das heißt, es muß 

gelten mit der Nukleonenzahl. Wendet man dieses Kriterium auf den Ausdruck oben 
an, so liegt für 7 > | ein Minimum von Evor, das heißt stabiles Gleichgewicht: 

1/(7-4/3) 

(2.48) 



3a(7 - 1) 



7 < I liefert instabiles Gleichgewicht. Der Fall 7 = | erfordert a = b, damit verschwindet 
die Energie. 
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Anwendung auf Weiße Zwerge 

Die relativistisclie kinetisclie Energie eines Teilchens ist gegeben durch 



^kin =Vp^ + M^-M . (2.49) 

Wir nehmen an, daß die Temperatur gegenüber der kinetischen Energie des Teilchens 
vernachlässigbar sei. Die Dichte der kinetischen Energie e und der Druck der Elektronen 
p in einem Weißen Zwerg sind gegeben durch (Rechnungen ab sofort in natürlichen 
Einheiten, d.h. h = c = 1) 

y (yPTmf-me) (2.50) 



(27r)3 





k 



F 



5 TT 



P = — TT / — fc^dfc 
3(27r)^y Jk^ + ml 
^ 

; /o 2 M/3 f 3vr2p '^B 

\mAr/i / ">-!^ 

Tie 

kp '■ Permi- Impuls, rzg : Anzahl der Baryonen pro Elektron . 
Zwei Extremfälle sind leicht zu diskutieren: kp -C mg und kp ^ me- Im ersten Fall 

P = — 77rTi4 = 7^ j ; e = -p . 

Das ergibt unter Verwendung von Gl. ()2.45|) eine Masse von 



M = 2.79 



Pc ) ml 



Pc : Dichte, wo der Fermi-Impuls gleich der Elektronmasse ist 

Im zweiten Fall kann die Elektronmasse überall vernachlässigt werden: 

Stt ,4 1 / 37r2 
P = T7T77rT?4 = TTTT P^! ^ = 3p 



12(27r)2 ^ 127r2 \mNP ) 
M = 5.87 — ^ Chandrasekhar-Grenze 



2.51) 

2.52) 

2.53) 
gilt: 

2.54) 



2.55) 
2.56) 

2.57) 
2.58) 
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In beiden Fällen haben wir es mit Polytropen zu tun: Im ersten Fall mit einem Adiabatenexponenten 
von 7 = 5/3, im zweiten Fall ist 7 = 4/3. Weiße Zwerge können unter diesen Voraussetzungen 
der Gravitation standhalten und sind bis zu der Grenzmasse in Gl. ()2.58j) stabil. 

Materie verhält sich aber etwas anders. Ab einer gewissen Größe des Fermi-Impulses 
wird Elektroneneinfang energetisch günstiger, so daß die Anzahl der Baryonen pro Elektron 
wächst. Dies wiederum reduziert die Grenzmasse. Man kann zeigen, daß an diesem Punkt 
Stabihtät in Instabihtät umschlägt (s. |Wei72] ). 

2.3 Endzustände der Sternentwicklung 

In diesem Abschnitt werden Endzustände der Sternentwicklung erläutert. Es gibt genau 
drei Endstadien der Sternentwicklung: 

• Weißer Zwerg 

• Neutronenstern 

• Schwarzes Loch 

Der Endzustand eines Sterns hängt von seiner Anfangsmasse ab und von den äußeren 
Lebensbedingungen. So können Materieströme von einem Begleiter des Sterns seine Masse 
erhöhen und dafür sorgen, daß die Entwicklung weitergeht. 

2.3.1 Weiße Zwerge 

Irgendwann ist der gesamte Brennstoff eines Sterns aufgebraucht. Der thermische Druck 
kann der Gravitation nicht weiter standhalten. Bei kleinen Massen verhindert die elektrostatische 
Abstoßung der Elektronen einen Kollaps. Ab einer Dichte von 80 g ■ cm~^ tritt wie beim 
Roten Riesen Entartung der Materie ein. Der Stern kühlt ab. Ist die Masse des Gores 
kleiner als die Ghandrasekhar-Grenze, und gibt es auch keine allzu schwere Hülle, so 
entwickelt sich der Stern zu einem Weißen Zwerg. Sein Radius beträgt einige tausend 
Kilometer. Die Elektronen im Inneren des Gores sind relativistisch entartet und können 
dadurch den Kollaps des Sterns verhindern. Die äußeren Schichten des Gores sind weniger 
dicht, und die Elektronen sind nicht entartet. Diese äußere Schicht wirkt als Isolation 
gegenüber dem Außenraum. Durch das Schrumpfen gewinnt der Stern Energie aus dem 
Gravitationsfeld und strahlt diese Energie in langen Zeiträumen ab, der Weiße Zwerg 
strahlt weiß-bläulich. 

Hat der Stern dagegen eine schwere Hülle, so kann der Kern sich so weit aufheizen, 
daß Helium zu Kohlenstoff verbrannt wird. 
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2.3.2 Neutronensterne 



Ist das Gore des Sterns schwerer als die Chandrasekhar-Grenze, wird der Stern instabil. 
Elektroneneinfang wird jetzt energetisch günstiger, wodurch sich die Zustandsgieichung 
der Materie ändert und ein neues Gleichgewicht erreicht wird. In Gl. fl2.58p kann für leichte 
Neutronensterne die Elektronenmasse rrie durch die Neutronenmasse m„, ersetzt werden, 
und rie wird gleich 1 gesetzt. Bei schwereren Neutronensternen bricht die Analogie zum 
Weißen Zwerg zusammen, weil die Materie eines Weißen Zwergs aus nichtrelativistischen 
Nukleonen besteht, während die Nukleonen eines Neutronensterns (von einer Sonnenmasse) 
eine kinetische Energie im Bereich ihrer Ruhmasse haben. Neutronen sind gegen ß- 
Zerfall stabil, wenn der Fermi-See der Elektronen im Gore mindestens bis zum maximalen 
Elektronenimpuls der beim /3- Zerfall entstehenden Elektronen aufgefüllt ist. Dieser maximale 
Impuls beträgt etwa 1.2 MeV ( |Wei72j ). Im Gore sind bei diesen hohen Dichten die 
Elektronen entartet. Ihr chemisches Potential beträgt am Anfang etwa 6 MeV ( |HHAL79j ) . 
Es bildet sich ein /9-Gleichgewicht aus, so daß am Ende das chemische Potential der 
Elektronen gleich der Differenz der chemischen Potentiale von Neutronen und Protonen 

ist: IJ,e = fln- (J'p- 

Im letzten Schritt wird bei massiven Sternen Silizium zu Eisen verbrannt. Dies geschieht 
innerhalb eines Tages. Wenn der thermische Druck wegfällt, schrumpft der Stern. Dabei 
werden die Elektronen in die Atomkerne gepreßt. Ab einer bestimmten Dichte wird 
Elektroneneinfang für den Stern energetisch günstiger. Dabei werden Neutrinos freigesetzt. 
Es werden aber nicht alle Protonen in Neutronen umgewandelt, sonst wäre der Stern 
instabil gegenüber /3-Zerfall. In dieser Umgebung können sogar Myonen und andere exotische 
Teilchen wegen des Pauli-Blocking stabil sein! Elektroneinfang reduziert den Permi- Impuls 
und somit den Druck und leitet einen Kollaps des Kerns ein, der von einer Explosion 
begleitet sein kann. Dabei wird die äußere Hülle weggesprengt. Dieses Ereignis ist eine Typ 
Il-Supernova. Übrig bleibt ein extrem dichter Neutronenstern, sofern der Stern genügend 
Materie losgeworden ist. Dieser wird von den nichtrelativistisch entarteten Neutronen am 
weiteren Kollaps gehindert. Die Rolle der Elektronen haben die Neutronen übernommen. 
Die Rechnung von Ghandrasekhar kann im wesentlichen übernommen werden, man muß 
nur Elektronen durch Neutronen ersetzen. Dies führt zu einer Grenzmasse in derselben 
Größenordnung. Der resultiertende Stern ist kleiner und dichter. 

Außen ist der Stern von einer festen Kruste aus Protonen in fester Phase umgeben. Die 
Elektronen sind relativistisch entartet und werden von den Atomkernen nicht beeinflußt. 
Im Inneren des Sterns befindet sich eine "Neutronenfiüssigkeit" . 
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2.3.3 Schwarze Löcher 

Ist der Stern zu massiv und schafft er es nicht, genug Materie los zu werden, kann der 
Entartungsdruck der Fermionen der Gravitation nichts mehr entgegensetzen. Der Stern 
kollabiert zu einem Schwarzen Loch. 

2.4 Sternkollaps und Neutrinokühlung 

In diesem Abschnitt wird die Explosion (Typ Il-Supernova) eines massiven 20 — 25Mq) 
gealterten Sterns genauer betrachtet (nach |BB85p . 

Beginn des Kollaps 

Wir beginnen mit dem Zeitpunkt, wo die Verbrennung von Silizium zu Eisen einsetzt. 
Dieser Prozeß dauert etwa einen Tag, während die anderen Prozesse in längeren Zeiträumen 
stattgefunden haben. Innerhalb dieses Tages bildet sich ein Gore aus Eisen, der die 
Ghandrasekhar- Grenze übersteigt. Dieses Gore kann nicht mehr durch den Elektronendruck 
stabilisiert werden, es implodiert zuerst, um dann zu explodieren. 

Wenn die Implosion einsetzt, wird das Gore komprimiert. Das führt zu einer Temperaturerhöhung, 
aber nicht unbedingt zu einem höheren Druck. Der Grund ist folgender: Der Druck hängt 
ab von der Anzahl der Teilchen und deren mittlerer Energie (ideales Gas: Druck ist 
proportional zur Teilchendichte und zur Temperatur). Im Gore befinden sich Elektronen 
und Eisenkerne. Da ein Eisenkern aus 56 Nukleonen aufgebaut ist und weniger Eisenkerne 
als Elektronen zur Verfügung stehen, ist der Beitrag der Elektronen zum Gesamtdruck 
wesentlich größer als der Beitrag der Eisenkerne. Durch das Aufheizen werden einige 
Eisenkerne in kleinere Bruchstücke zerlegt, dieser Vorgang benötigt soviel Energie, wie 
bei deren Produktion freigesetzt wurde. Diese Energie wird den Elektronen entzogen, 
wodurch deren Druck sinkt. Durch das Aufspalten dieser Eisenkerne steigt der Druck 
der Nukleonen, aber bei weitem nicht so stark, daß der Druckabfall durch die kühleren 
Elektronen kompensiert werden kann. Die Folge dieses Druckabfalls ist eine Beschleunigung 
des Kollapses. 

Die Entropie 

An dieser Stelle ist eine Bemerkung zur Entropie angebracht. Der Kollaps ist kein chaotischer 
Prozeß, wie es auf den ersten Blick scheint. In einer Wasserstoffwolke ist die Entropie hoch, 
die ja ein Maß für die Unordnung ist. In dem Gore eines kollabierenden Sterns sind aber 
56 Nukleonen zu einem Eisenkern geordnet, folglich ist die Entropie pro Nukleon klein. 
Photonen und Neutrinos haben im Laufe der Zeit die fehlende Entropie abgeführt. Da 
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die Zeitskala des Kollaps in der Größenordnung von Millisekunden liegt, Starke Prozesse 
aber eine Zeitskala von 10~^°s haben, werden Störungen sofort beseitigt, weil der Kollaps 
im Vergleich zu der Zeitskala starker Prozesse "ewig" dauert. 

Bedingt durch die hohe Dichte und der damit einhergehenden Entartung der Elektronen 
tritt vermehrt Elektroneneinfang ein. Es enstehen dabei Neutrinos, die bei ihrer Emission 
Energie und Entropie abführen. Der Gesamteinfluß auf die Entropie ist unsicher. 

Die verminderte Elektronenzahl reduziert weiter den Druck und beschleunigt so den 
Kollaps. Da das Verhältnis der Anzahl von Elektronen zu Baryonen kleiner wird, wird 
auch die Chandrasekhar- Grenze kleiner, welche ja von diesem Verhältnis abhängt. 

Neutrinoopazität 

Bei einer Dichte von 4 ■ 10^^ g • cm~^ tritt ein Effekt auf, der sonst nur unmittelbar 
nach dem Urknall zu beobachten gewesen wäre: Die Sternmaterie wird für Neutrinos 
undurchsichtig. Die mittlere Weglänge A eines 30 MeV- Neutrinos beträgt im Gore etwa 
300 cm. Die Diffusionszeit (Zeit, die ein Teilchen braucht, um das Gore zu verlassen) eines 
solchen Neutrinos beträgt damit etwa Is (nach |Raf96j ). Die Neutrinos können das Gore 
nicht verlassen und es nicht kühlen. Hier ändert sich auch die Rolle der Ghandrasekhar- 
Grenze. Vorher war das die größte Masse, die ein entartetes Elektronengas am Einsturz 
hindern kann. Jetzt ist das die Grenzmasse, die als Ganzes kollabieren kann. Gebiete 
innerhalb dieses Teils des Gores können durch Schall- und Druckwellen erreicht werden 
und somit Dichteschwankungen ausgleichen. 

Aufbau einer Schockfront 

Der Einsturz des Gores geht weiter. Irgendwann hat die Materie im Zentrum Kerndichte 
und darüber (pKern2.7 ■ 10^^ g ■ cm) erreicht, weil von oben immer neues Material einfällt. 
Kernmaterie ist zwar schwer komprimierbar, aber es ist möglich. Die inneren Schichten 
stehen dadurch unter einer großen Spannung wie eine zusammengedrückte Feder. Die 
Implosionsfront kommt mit einer gewaltigen Erschütterung zum Stillstand und wird zu 
einer Explosionswelle, die durch die Entspannung des Gores ausgelöst wird. Die Geschwindigkeit, 
mit der eine Goreschale einstürzt, ist proportional zu ihrer Entfernung zum Sternzentrum. 
Man sagt deshalb, der Einsturz sei homolog. Die Dichte nimmt nach außen hin ab, die 
Schallgeschwindigkeit ebenfalls, da die Kompressibilität mit kleiner werdender Dichte 
zunimmt. Im Gore gibt es eine Grenze, wo die Einfallgeschwindigkeit gleich der Schallgeschwindigkeit 
ist. Diese markiert die Grenze des homologen Gores und wird sonischer Punkt oder 
Schallpunkt genannt. Eine Störung innerhalb des homologen Gores hat keinen Einfluß 
auf Gebiete außerhalb dieses Punktes, denn eine Schallwelle, die am Schallpunkt durch 
die einfallende Materie nach außen läuft, bleibt in bezug auf das Sternzentrum in Ruhe, 
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da dieses Coreelement mit Schallgeschwindigkeit einstürzt. Demzufolge kann eine Störung 
im Inneren den Außenbereich nicht erreichen. 

Wenn die Implosionsfront zum Stillstand gekommen ist, entstehen durch die Erschütterung 
Schallwellen, die nach außen laufen, mit größer werdender Entfernung zum Zentrum 
langsamer werden und am Schallpunkt schließlich stehen bleiben. Auf die Materie im 
Zentrum, die ja Kerndichte erreicht hat, fällt immer mehr Material und erzeugt weitere 
Schallwellen, die sich am Schallpunkt sammeln. Dort baut sich eine Druckfront auf. Diese 
bremst auffallende Materie, wodurch ein Sprung in der Schallgeschwindigkeit entsteht. 
Diese Unstetigkeit verursacht eine Schockwelle. 

Nachdem das Gore Kerndichte erreicht hat und durch weiteres einfallendes Material 
weiter komprimiert worden ist, entspannt sich das Gore und sendet weitere Schallwellen 
aus, die sich mit der Schockfront am Schallpunkt vereinen. 

Schockwellen 

An dieser Stelle soll der Unterschied zwischen einer Schall- und einer Schockwelle geklärt 
werden. Schallwellen durchlaufen ein Medium mit Schallgeschwindigkeit und verändern 
die Eigenschaften des Mediums nicht. Es wird vor allen Dingen nicht zerstört. Eine 
Schockwelle pflanzt sich mit Uberschallgeschwindigkeit fort und ändert dabei Druck, 
Dichte, Entropie des durchlaufenen Mediums. 

Neutrinokühlung 

Die bei der Supernova SN 1987A beobachteten Neutrinos wurden in einer Zeitspanne von 
etwa 10 s emittiert. Nach etwa einer Millisekunde werden die Neutrinos freigesetzt, die 
beim Elektroneneinfang entstehen. Dies ist der prompte Neutrinoausbruch. Es bleiben 
aber noch genügend Leptonen (Neutrinos und Elektronen) im Inneren des Gore gefangen, 
wo sich die Schockwelle gebildet hat. Der äußere Teil des Gores fällt in der nächsten 
Sekunde in sich zusammen. Seine Bindungsenergie wird in Form von Neutrinos emittiert. 
Mittlerweile konnte die Schockwelle weiter nach außen laufen und die äußeren Schichten 
des Sterns wegreißen. Das Gore (Radius etwa 30 km) kann als Stern im Stern betrachtet 
werden. Er kontrahiert weiter und wird durch die Emission von Neutrinos aller Generationen 
gekühlt. Durch die Emission von Elektronneutrinos gehen die Leptonen bis auf einen 
kleinen Rest verloren. Nach 5—10 s ist diese Phase abgeschlossen (sogenanntes Helmholtzkühlen, 
s. [Raf96 ). Danach setzt die späte Kühlungsphase des Neutronensterns ein. 

Der wesentliche Punkt ist die Dauer der Neutrinokühlung von ^ 10 s, die im Vergleich 
zur Zeitdauer des Kollaps sehr lang ist. Die Neutrinos können aus dem Gore nicht frei 
herausströmen, denn ihre mittlere Weglänge ist kleiner als der Goreradius. Sie müssen 
aus dem Gore langsam herausdiffundieren. Dieses Neutrinosignal ist 1987 nachgewiesen 
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worden (IMB, Kamiokande) . Ein zusätzlicher Kühlmechanismus darf deshalb nicht zu 
effektiv sein. Für die mittlere Weglänge gilt nämhch X^, — l/{na) oc l/(nT^), wobei n 
die Teilchendichte, er den totalen Wirkungsquerschnitt und T die Temperatur bezeichnen. 
Wäre ein Kühlmechanismus effektiver, würde die Coretemperatur rasch sinken und die 
mittlere freie Weglänge der Neutrinos größer werden. Ein kürzeres Neutrinosignal wäre 
die Folge. Das steht aber im Widerspruch zu den Beobachtungen. 

Weitere Verlauf 

Das Geschehen danach ist unklar. Im einfachsten Szenario läuft die Schockwelle durch 
alle Schichten des Sterns in Richtung Oberfläche: Der Stern explodiert, man kann eine 
Supernova beobachten. Der größte Teil der Hülle wird weggesprengt, übrig bleibt ein 
kleiner, extrem dichter Reststern, der zu einem Neutronenstern kondensiert. Die Schockwelle 
setzt auf ihrem Weg nach außen die gefangenen Neutrinos frei, die den Stern dabei kühlen. 

Computersimulationen zeigen aber ein anderes Bild. Die Schockwelle läuft nach 100 — 
200 km aus, weil sie auf ihrem Weg die Atomkerne in kleinere Bestandteile (z. B. Protonen) 
zerlegt und dabei an Energie verliert. Druck und Temperatur werden stark vermindert. Die 
Protonen fangen Elektronen ein, die dabei entstehenden Neutrinos entkommen und führen 
Energie ab. Wegen der vielen Zwischenfälle, die der Schockwelle in dieser Region passieren 
können, wird dieser Teil des Sterns "Minenfeld" genannt. Es gibt kein befriedigendes 
Modell, das die Explosion zutreffend beschreibt. 

Das soeben beschriebene Szenario ist ein mögliches unter vielen. Simulationen haben 
gezeigt, daß der Mechanismus stark von der Masse und der Kompressibilität der Kernmaterie 
abhängt. Das letztgenannte Modell trifft wohl auf Sterne mit einer Anfangsmasse von 
12 - 18 Mq zu. 

Schicksal des Reststerns 

Zum Schluß noch ein paar Bemerkungen zu dem Schicksal des Gores. Die Gravitation kann 
den Elektronendruck überwinden. Sogar Kernmaterie kann den weiteren Einsturz nicht 
immer verhindern. Die Grenzmasse für einen erkalteten Neutronenstern liegt bei etwa 
1.8 Mq. Die Argumente sind dieselben wie bei der Chandrasekhar-Grenze. Die nicht- 
relativistisch entarteten Neutronen können den weiteren Kollaps bis zu dieser Grenze 
verhindern. Größere Massen kollabieren wahrscheinlich zu einem Schwarzen Loch. Unmittelbar 
nach der Explosion ist der Reststern in seinem Inneren etwa 10^^ K heiß. Diese Temperatur 
reicht aus, um größere Massen als die 1.8 zu stabilisieren. Die Kühlung erfolgt durch 
Neutrinos. Diese führen hundert Mal mehr Energie als die Explosion ab, in Zahlen ^ 
lO^^erg, was etwa einem Zehntel der Gesamtmasse des Sterns entspricht. 
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Auf der Erde können diese Modelle durch Nachweis der Neutrinos und der Häufigkeit 
der chemischen Elemente überprüft werden. 



Kapitel 3 

Neutralinokühlung der Supernova 




3.1 Die Emissivität des Elektron-Positron-Kanals 

In diesem Kapitel werden die Produktionsraten für Neutralinos durch die verschiedenen 
Mechanismen berechnet. 

3.1.1 e''"e~- Annihilation und Selektronaustausch 

Wir nehmen an, daß folgende e+e^-Prozesse zur Kühlung der Supernova beitragen können: 



Der erste Prozeß ist nur dann möglich, wenn zum leichtesten Ncutralino Higgsino-Anteile 
beigemischt sind, weil mir die Higgsinos an das koppeln. Die Feynman-Graphen für 
diese fünf Prozesse sind in erster Ordnung: 




^ ^ XiXi 



(3.1) 

(3.2) 




e ,P2 



Xi,P4 e ,p2 




Xl,P4: ^- 



Xl,P4. 




e^,Pi 



e^,Pi 



Xl,P3 ^^Pi 



Xl,P3 



Xi,P3 



29 



30 



KAPITEL 3. NEUTRALINOKÜHLUNG DER SUPERNOVA 



Mit den Feynmanregeln des Standardmodells (ich beziehe mich auf |HM84j ) und des 
MSSM aus |HK85j und |(7H86j kann man den totalen Wirkungsquerschnitt berechnen. 
Bartl et al. haben in |BFM86j diesen Wirkungsquerschnitt in allgemeinerer Form für zwei 
verschiedene Neutralinos als Endzustand berechnet. Da die anderen Neutralinos zu große 
Massen haben als daß sie thermisch angeregt werden können, beschränke ich mich auf 
den Sonderfall, daß zwei leichteste Neutralinos auslaufen. Für die ausführliche Rechnung 
verweise ich auf den Anhang, hier nenne ich nur die wesentlichen Ergebnisse. 

Es ergeben sich folgende Amplituden: 

-JL^Dzo{s)v{pi)j''{clPl + CRPn)u{p2Mpsh^{LiPL + RiPr)v{p^.3) 

-g%\^D,{t)v{pi)PRvips)ü{p^)PLu{pi) (3.4) 
-g%\''D,{u)vipi)Pnvipi)üip3)PLuip2) ■ (-1) (3.5) 



-iMi 

-iM2 

-iMs 
-iM4 
-iMs 



-g^\ri\^ De(t)v{pi)PLv{p3)u{p^)PRu{p2j 

-g'^\ri\'^De{u)v{pi)PLv{pi)u{p^)PRu{p2) ■ 



(3.7) 



Summation über die Ausgangsspins und Mittelung über die Eingangsspins liefert den 
differentiellen Wirkungsquerschnitt, ausgewertet im Schwerpunktssystem: 



da 



9" 



'1- 



^ ( (t - M'f + {u- My - 2Mh] X 



■2 ^1 Li{cL\li\^-CR\n\') \h\'+\n\' 



die auftretenden Konstanten bedeuten: 



M|o cos4 öw ^ Ml^Mj cos2 Ow ^ 4MJ 



(3.8) 
(3.9) 



1 1 
cl = T^icv + Ca); cr = -(cy - Ca); cy = T3- 2ee sin^ 9w; ca = T3 



2 

^ Ni,N*,,) sin(2/?) 



Li = liN,,N*,-Ni,N*,)cosm-liN,,N*,-Ni, 
Ri = —Li 
h = V2 



i sin Qw^ii 



1 



'1 



cos Qw \ 2 



- — Ce Sin^ Q^JY ) A''i2 



ri 



= 72 



iVi 



• a AT* ^ 
cos 



DAS) = D,{x) = j^^ 



(3.10) 

(3.11) 
(3.12) 

(3.13) 
(3.14) 
(3.15) 
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Nach Integration über die Winkel erhält man den totalen Wirkunsquerschnitt: 



a 



1-na s 
3 sin^ e 



w 



L\^c\L\ , Li(cz.|/ip-c^|ri|2) |Zi|4+|ri 



M|o cos4 dw 



+ 



M|oM| cos2 



+ 



ml 



(3.16) 



Weitere Details und Erklärungen zu dieser Rechnung befinden sich im Anhang. 



3.1.2 Definition der Emissivität 

Die Emissivität e gibt den Energieverlust pro Zeit und Volumen an. Sie kann analog 
zum Wirkungsquerschnitt aus dem Betragsquadrat des Matrixelements unter Verwendung 
Fermis Goldener Regel berechnet werden. Sie ist für einen 2 2-Prozeß gegeben durch: 



(2vr)2U (27r)4E2 {^-aflE, {l-aflE,^^' + ^4)/i/2(l - /3)(1 - A) 
^\vx^V2-Vz-vi)\Mf (3.17) 
fi = ^p,^^,yrp ^ ^ (alle (Anti-) Teilchen sind Fermionen). (3.18) 

Aus der Formel ist der Zusammenhang zwischen Emissivität und Wirkungsquerschnitt 
ersichtlich: 

e = j^jdpsdp,{E3 + E^)M2{l-h){l-h)vu0ia . (3.19) 

An dieser Stelle sei noch darauf hingewiesen, daß die beiden einlaufenden Teilchen im 
Sternsystem nicht kollinear sind. Deswegen muß für |vi — V2I die M0ller-Geschwindigkeit 
angesetzt werden (siehe |LL71j ). die als Funktion der Teilchengeschwindigkeiten folgende 
Form hat: 



V(V1-V2)2-(V1 X V2)2 . (3.20) 



Um die Emissivität ausrechnen zu können, muß man die Temperatur und das chemische 
Potential kennen. Wenn die Temperaturen hoch und damit die Entartung niedrig ist, so 
sind /a und Null. Die Faktoren 1 — /s, 1 — /4 geben die Phasenraumunterdrückung an, 
da Fermionen nur einen freien Zustand pro Teilchen besetzen dürfen, fi und /2 geben die 
Teilchendichte an. 
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3.1.3 Emissivität des Elektron-Positron-Prozesses 



Der totale Wirkungsquerschnitt des Prozesses e+e" — > wurde schon ausgerechnet. 
Das Ergebnis kann jetzt für die Berechnung der Emissivität verwendet werden. Die Impulse 
der einlaufenden Teilchen (Elektron und Positron) werden wie folgt parametrisiert: 



Pl 



(j9i,0,0,pi), = (p2, 0,P2 sin 0,^2 cos I 

^(Pl,P2) 



(3.21) 
(3.22) 



Einige Winkelintegrationen in Formel ()3.19j) können schon ausgeführt werden. Man erhält 
als Resultat, wenn man Ei^2 = Pi,2 (masselose Elektronen bzw. Positronen) und E1 + E2 = 
E3 + E4 (Energieerhaltung) benutzt (siehe [KacOnj ) : 



dpidp2 dcos 9- 



p\ 



-{Pi +P2)vM0\(y 



(3.23) 



Unter Verwendung von Gl. ()3. 2111 und v = p/i? erhält man für die Möller- Geschwindigkeit: 



■^1^01 



|Vl — V2IM0I 



Pl 


P2 




Pl 
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'1 — cos I 



(3.24) 



3.1.4 Die abgestrahlte Gesamt energie 

Integration der Emissivität über das Sternvolumen und die Abstrahldauer liefert die 
Leistung beziehungsweise die Energie dieses Prozesses: 



E 
At 



dt j dVe{M^,p,T,7^)); 



At AV 

Abstrahldauer AV : Sternvolumen 



(3.25) 
(3.26) 



Weil in |BL86j Temperatur und Entartung (siehe die entsprechenden Abb. aus |BL86j ) 
nicht in Abhängigkeit des Radius, sondern der eingeschlossenen Baryonenmasse angegeben 
sind, muß dV^ auf dm umgerechnet werden. Das Dichteprofil aus jBL86j zeigt, daß die 
Dichte des kollabierenden Gores nach 250 ms etwa konstant ist. Damit ist die Umrechnung 
einfach: dV = d— = —. Für die Abstrahldauer nehme ich Zeitspannen zwischen 0.005 s 
und 2 s an, weil die Zeitspannen, in denen der Kollaps stattfindet, diesen Größenordnungen 
entspricht. Für die anderen Parameter fixiere ich sie bei 1 s. 



3.1. DIE EMISSIVITÄT DES ELEKTRON-POSITRON-KANALS 



33 



3.1.5 Temperatur und chemisches Potential 

Ein Ziel meiner Arbeit ist es, für die Temperatur und das chemische Potential eine 
realistische Verteilung anzunehmen. Burrows und Lattimer haben in |iBL86j verschiedene 
Parameter einer Supernova untersucht, unter anderem haben sie Profile für die Temperatur, 
chemische Potential und Dichte zu verschiedenen Zeiten berechnet. Diese Profile verwende 
ich, um die Emissivität in Abhängigkeit des Sternradius zu berechnen. Ich benutze die 
Profile zur Zeit t = 0.5 s, weil diese in der Mitte der ersten Sekunde liegen. Aus den 
Funktionsgraphen werden etwa ein halbes Dutzend charakteristische Punkte (z. B. Extrema) 
abgelesen, und diese dann durch kubische Splines miteinander verbunden. Um ein gutes 
glattes Resultat zu erzielen, reichen die wenigen Punkte aus. Ein Problem bei diesen 
Simulationen ist die Unkenntnis der Anfangsdaten. Die Simulation von Burrows und 
Lattimer beginnt mit dem Kollaps des Sterns. Die Simulationen können aber noch nicht 
den Prozeß des Kollabierens erfolgreich berechnen. Die Schockwelle bleibt im Gore stecken. 
Deshalb sind die Anfangsdaten nur unzureichend bekannt. Nach einer gewissen Zeit 
s) steht das Supernovasystem aber im Diffusionsgleichgewicht, so daß die genaue 
Kenntnis der Anfangsdaten für die weitere Entwicklung nicht nötig ist. Man erkennt 
aus den Graphen, daß bei r ^ Mq die Temperatur hoch und die Entartung gering 
ist. Die numerische Auswertung der Emissivität zeigt, daß die meiste Energie aus einer 
Kugelschale mit diesem Radius herrührt. 

Der Radius ist als eingeschlossene Baryonenmasse angegeben. Uber die Goredichte 
kann dann die Massenabhängigkeit in eine Radiusabhängigkeit umgewandelt werden. Aus 
der Größe des Sterns R = 13km und seiner Masse M = 1.4 Mq ergibt sich p zu 

mm 14 ■? / N 

P=^ = i^^S-lO^^gcm^ ; (3.27) 

das Gore hat also Kerndichte erreicht. 

Alle zur Berechnung der Emissivität notwendigen Größen sind jetzt beisammen. Die 
drei verbleibenden Integrationen in Gl. ()3.23j) werden numerisch ausgeführt. Dazu verwende 
ich die verallgemeinerte Gaußsche Integrationsmethode, bei der der Integrand f{x) an den 
Nullstellen orthogonaler Polynome ausgewertet wird und diese Werte gewichtet addiert 
werden. Die zu verwendenden Polynome hängen von der Gestalt des Integranden ab: 

b 

/n 
k{x)f{x)dx ^ Y^afixi), (3.28) 

k{x) : nichtnegative, meßbare Gewichtsfunktion, (3.29) 
Cj : Gewichtsfaktoren . (3.30) 
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Die Intervallgrenzen können auch uneigentlich sein. Die Xj sind Nullstellen von Polynomen, 
diese Polynome bilden bezüglich des Skalarprodukts 

b 

U.9) = j k{x)f{x)g{x)dx (3.31) 

a 

eine Orthonormalbasis. Die genaue Theorie dazu findet man z. B. in |Sto94j . 

Für dieses Problem geeignete Polynome sind für die Impulsintegrationen die Laguerre- 
Polynome, da die Fermi-Dirac- Verteilungen sich wie ö{e~^) verhalten, für die Winkelintegration 
wähle ich die Legendre-Polynome. Die Nullstellen Xj und Gewichtsfaktoren Cj findet man 
in |SS66j oder in |AS70j tabelliert. Dieses Verfahren hat den Vorteil, daß man mit wenigen 
Schritten (und damit in kürzerer Zeit gegenüber Keplerscher Faßregel oder ähnlichem) ein 
genaues Endergebnis erhält, nämlich die Emissivität in Abhängigkeit von der Masse und 
der Kopplung, das heißt, den Mischungsparametern und der Squarkmasse. 



3.1.6 Das Resultat 

In Abbildung 13.11 ist die Leistung in Abhängigkeit der Neutralinomasse exemplarisch für 
ein Bino mit einer Selektronmasse von 100 GeV dargestellt. 

Zur besseren Handhabung der Funktion habe ich sie durch eine geeignete analytische 
Funktion approximiert. Folgender Ausdruck ist eine mögliche Näherung: 



P{M,) = 6.47 ■ 10^ j exp -0.0115 ^ - 0.000053 j j 10^^ erg/^3.32) 



Die Näherung ist als gestrichelte Funktion ebenfalls in Abb. l3.ll zu sehen. Diese Näherungsformel 
hat bis = 370 MeV einen maximalen Fehler von 20%, darüber hinaus ist der Fehler 
so groß, daß die Formel nicht mehr brauchbar ist. 



3.1.7 Frühere Resultate 

EUis et al. haben in [EOSSSg] die Möglichkeit von Photinos mit einer Masse von etwa 
100 eV diskutiert. Squark- und Selektronmassen in Bereich der VT-Masse zerstören eventuell 
den durch Neutrinobeobachtungen erlaubten Bereich. Neutralinos sind dann genauso wie 
Neutrinos gefangen, während Sfermionmassen von 2.5 TeV die Photinos frei herausströmen 
lassen. 
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Leistung in Abhängigkeit der Neutralinomasse 

(Bino-ähnlich, M- = 100 GeV) 




2 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

100 200 300 400 

(MeV) 



Abbildung 3.1: Leistung in Abhängigkeit von der Masse für ein Bino-Neutralino, Mg = 100 GeV 

Kachelrieß hat in [KacOnj die Möglichkeit eines 34 MeV-Neutralinos untersucht und 
gefunden, daß die relative Emissivität (Emissivität normiert auf die zentrale Emissivität) 
nur schwach unterdrückt ist. 

Dicus et al. haben in |l)M^r98| eine unter Massenschranke für ein superleichtes Gravitino 
mittels Supernova-Kühlung ausgerechnet. Sie haben als maximale Leistung lO^^erg/s 
angenommen und damit die Gravitino-Emission nach oben begrenzt. 

3.1.8 Diskussion der Parameter 

Abschätzung der abgestrahlten Energie 

Um eine untere Schranke für die Neutralinos zu gewinnen, braucht man eine obere Schranke 
für die maximal zulässige abführbare Energie, ohne daß das Neutrinosignal gestört wird. 
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Einfache Überlegungen führen zu folgender Formel für die Bindungsenergie, die beim 
Einsturz einer kugelsymmetrischen Massenverteilung frei wird (aus ^R.af96] ) : 



Nimmt man einen Endradius des Gores von 13 km wie Burrows und Lattimer in 
|BL86j an, so erhält man als freiwerdende Bindungsenergie 3 • 10^^ erg. Diese Energie 
wird in Form von Photonen und Neutrinos abgegeben. Obwohl eine Supernova eine 
ganze Galaxie überstrahlen kann, ist die abgestrahlte Energie in Form von Photonen 
gegenüber den Neutrinos vernachlässigbar. Es werden etwa 2.4 — 3.0-10^^ erg als Neutrinos 
abgestrahlt, das entspricht etwa 50000 Erdmassen oder 0.15 Sonnenmassen f |Bur90j )! Da 
es sechs Neutrinoarten gibt, führt jede im Mittel 4 — 5 ■ 10^^ erg ab. Damit kann ich 
die maximale Energie, die die Neutralinos abführen können, zn E = 10^^ erg mit einer 
Unsicherheit von der Größenordnung 0{1) festlegen. In Abb. wird die Abhängigkeit 
der unteren Massenschranke von der erlaubten abgestrahlten Energie für verschiedene 
Selektronenmassen dargestellt. 

Für größere Energieschranken fällt die untere Massenschranke für das Neutralino 
bei konstanter Selektronenmasse schwach monoton, wie man das auch erwartet. Die 
Selektronenmasse hat einen erheblichen Einfluß. So liegt die untere Massengrenze bei 
einem leichten Selektron (100 GeV) bei 350 . . . 400 MeV, bei einem schweren Selektron 
(1000 GeV) ist das Neutralino leichter: 100 .. . 200 MeV. 

Die Zeitskala 

Weiterhin muß man abschätzen, in welcher Zeit die Energie abgeführt wird. Der Kollaps 
flndet auf einer Zeitskala von Bruchteilen einer Sekunde bis eine Sekunde statt. Innerhalb 
einer Sekunde werden 10^^ erg abgeführt, damit hätte ein Bino (siehe Abb. 13. Ij) eine 
untere Massenschranke von etwa 360 MeV. Vergrößert man die Abstrahlzeit (die maximal 
abgestrahlte Energie wird konstant bei 10^^ erg gehalten) für verschiedene Selektronenmassen, 
so wird die untere Massenschranke für das Bino kleiner, wie man aus Abb. 13.31 ablesen 
kann. Bei einem Selektron mit einer Masse von 1 TeV ist hier keine Aussage möglich. Für 
sehr kleine Zeiten hängt die Massenschranke empflndlich von der gewählten Abstrahlzeit 
ab, weil in diesem Bereich die Kurven steil verlaufen. 

Die Mischungsparameter 

Betrachtet man die Abhängigkeit von den Mischungsparametern, so muß man im M — 
M' — IX — tan /3-Parameterraum nach Gebieten suchen, die folgende Bedingungen erfüllen: 




(3.33) 
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Abhängigkeit der Massenschranke von der Energieschranke 
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Abbildung 3.2: Abhängigkeit der Massenschranke von der maximal zulässigen abgestrahlten 
Energie 
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Massenschranke für ein Bino 

in Abhängigkeit von der Strahlungsdauer 
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Abbildung 3.3: Abhängigkeit der Massenschranke von der Dauer der Energieabstrahlung 
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10 MeV < M^o < 400 MeV 



My^± > 120 GeV (siehe einführendes Kapitel) 



13 + Nf^ < 0.5 = (0.08 



|l/4 



Diese Bedingung folgt aus dem Fehler der unsichtbaren Zerfallsbreite des Z'^-Bosons, 
siehe [CDRSOO, 



Eine ähnliche Suche im Parameterraum wurde in |(yDRSOÖ] durchgeführt, allerdings mit 
einer engeren Massenbreite und einer höheren Charginomasse. Für /x < GeV gibt es 
keinen zulässigen Eigenwert im Parameterraum. Den Parameter M2 muß man in sehr 
kleinen Schritten erhöhen, um das enge Tal zulässiger Werte nicht zu verfehlen. Da der 
kleinste Eigenwert von Interesse ist und "kleinster Eigenwert" einer Matrix, die von einer 
Reihe von Parametern abhängt, nicht unbedingt eine differenzierbare Funktion ist, kann 
der zulässige Bereich etwas sprunghaft verlaufen, wie die Abb. üOl zeigt. 



Bereiche mit Neutralinomasse 10 MeV < M < 400 GeV 

X 

|i = 300 GeV 



400 



tan(ß) = l 

tan(ß) = 8 

-- X^= 120 GeV, tan(ß) = 1 
120 GeV, tan(ß) = 8 




M' (GeV) 



Abbildung 3.4: Bereiche mit Neutralino zwischen 10 und 400 MeV, oberhalb der gestrichelten 
Linie ist die Masse des Chargino größer als 120 GeV 
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Für eine Selektronmasse in der Nähe der Z^-Masse ändert sich die Gesamtkopplung bei 
einer Variation der Mischungsparameter wenig, das heißt, man kann von der Massenschranke 
nicht auf möghche Mischungsparameter schheßen. In diesem Szenario ist der Z°-Beitrag 
äußerst klein. Das ändert sich, wenn die Selektronmasse hoch wird. Dann wird der Selektron- 
anteil klein gegenüber dem Anteil. Für zwei weitere mögliche Neutralinomischungen 
sind auf den folgenden Abbildungen die Massenschranken dargestellt. Die Energieschranke 
b( 



Neutralinoschranke 



> 



c 

C 



(Nj =0.9825, Nj2=-0.0847, Nj3=0.117, Nj4=-0.117, ^=300, tan(ß)=l) 



400 



300 




Zeit (s) 



Abbildung 3.5: Szenario 1 

Der Einfluß der Mischungsparameter A^i3 und A^i4 ist bei hohen Selektronmassen 
deutlich zu erkennen. Bei Szenario 1 sind die letztgenannten Parameter gleich groß: 
= Nu = 0.117, ^/NfPnWi = 0.165, im zweiten Fall ist A^ig = 0.221, Ar^4 = 
0.027, ^/Nfs + Nf^ = 0.223 (Dies ist der größte Wert, den ich für ^/Nf^ + Nf^ gefunden 
habe.). Im ersten Fall tritt partiell Auslöschung der Parameter auf, im zweiten Fall 
dagegen nicht. Außerdem spielt ein weiterer Effekt eine Rolle. Die Kopplungskonstante 
ist eine Summe aus Selektron- und Interferenzterm. Die Z°-Masse ist bekannt, die 
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Neutralinomassenschranke 

(Nj =0.975, Nj2=-0.024, Nj3=0.221, Nj4=-0.027, )i = 200 GeV, tan(ß)=8) 
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300 
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Abbildung 3.6: Szenario 2 



Selektronmasse nicht. Bei hohen Selektronmassen (TeV-Bereich) spielt der Selektronterm 
gegenüber dem Z°-Term keine Rolle mehr, da die Massen mit der inversen vierten Potenz 
eingehen, die Mischungsparameter jedoch nur quadratisch. Man erwartet also, daß die 
untere Massenschranke nicht mehr von der Selektronmasse abhängt, wenn das Selektron 
entsprechend schwer ist. Der Effekt der Mischungsparameter ist in den Abbildungen 13.51 
und l3.6l bei einer Selektronenmasse von 500 GeV deutlich zu erkennen, der andere dagegen 
nicht. Zwischen Szenario 1 (Abb. 13.5p und einem reinen Bino (Abb. 13. !( mit t = Is) ist 
kein Unterschied zu erkennen, da der Effekt der Mischungsparameter sich herauskürzt. 



Selektronenmasse 

Bisher ist die Selektronenmasse als dritter Freiheitsgrad bei den bisher diskutierten Parametern 
betrachtet worden. Nun werden diese fixiert (Energieschranke 10^^ erg, Bino) und die 
untere Massenschranke nur in Abhängigkeit der Selektronenmasse diskutiert, siehe Abb. 
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13.71 Wie zu erwarten war, verkleinert sich die untere Massenschranke bei größer werdender 
Selektronenmasse. Für noch größere Selektronenmassen ist meine Rechnung vermuthch 
nicht mehr brauchbar. Die untere Massenschranke variiert zwischen 38 MeV (Mg = 1400 GeV) 
lind 3a^MfiV (M. = IOOQpVV 

Abhängigkeit der unteren Massenschranke von der Selektronenmasse 

400 I 1 1 1 1 1 1 1 1 1 1 1 




Q I I I I I I I I I I I I 

500 1000 

Selektronenmasse (GeV) 

Abbildung 3.7: Abhängigkeit der unteren Massenschranke des Bino von der Selektronenmasse 
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3.2 Die Nukleon-Nukleon-Bremsstrahlung 

3.2.1 Die Mechanismen der Nukleon-Nukleon-Bremsstrahlung 

Es gibt einen weiteren Mechanismus, der Neutralinos erzeugen kann: Nukleon-Nukleon- 
Bremsstrahlung. Um die möglichen Reaktionen zu erhalten, nimmt man die Diagramme 
der e"'"e~-Prozesse und baut sie in jedes Beinchen der Nukleon-Nukleon-Streuung ein 
(Natürlich ist Elektron durch Quark und Selektron durch Squark zu ersetzen.). Die abgestrahlte 
Energie in Form von Neutralinos wird als weich betrachtet, obwohl dies für hohe Massen 
nicht zutrifft. Aber für eine Größenordnungsabschätzung ist diese Methode trotzdem 
angemessen. Die Diagramme für die Kopplung der Neutralinos an die Quarks sind in 
Tabelle ITT] aufgelistet. Später wird von den Quarks zu Nukleonen übergegangen. 





Tabelle 3.1: Nukleon-Nukleon-Bremsstrahlungsprozesse 



Die Feynmanregeln sind im Anhang 1X1 zusammengestellt. Die Matrixelemente haben 
die Form Pr ® bzw. ® Pr, falls ein Squark ausgetauscht wird, und {clPl + 
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crPr) (S> {LiPl + RiPr), falls ein Z°-Boson abgestrahlt wird. Weil die Sterntemperatur 
im Vergleich zur und zur Squarkmasse klein ist (30 MeV gegenüber fti 100 GeV), 
werden die Impulse in den Propagatoren vernachlässigt. Dadurch entsteht eine effektive 
Vierpunktwechselwirkung. Anschließend werden die Ströme durch eine Fierztransformation 
umgeordnet, so daß hadronischer Strom und Neutralinostrom miteinander wechselwirken. 
Unter Verwendung der S'f/(6)-Wellenfunktionen werden die Quarks in die Neutronen 
eingebunden. Dann wird durch den Neutralen Strom normiert, um nichtperturbative 
Effekte der Starken Wechselwirkung zu berücksichtigen. Dieser Prozeß ist in Tab. 13.21 
graphisch dargestellt. Die Neutralinos können an allen vier äußeren Beinchen des Graphen 
abgestrahlt werden. 



3.2.2 Der Neutralinotensor 

Näherungen 

Um die Bremsstrahlungsrate zu berechnen, nimmt man die Prozesse ohne Strahlung und 
fügt nacheinander an jedes Beinchen die Strahlung an. Die Nukleonen werden nichtrelativistisch 
behandelt, die Neutralinos dagegen relativistisch. Durch die nichtrelativistische Näherung 
für die Nukleonen vereinfachen sich die möglichen Kopplungsstrukturen, es gibt in diesem 
Fall nur zwei: den Einheitsoperator 1 und den Spinoperator er (mehr dazu siehe Anhang 
IB|1 . Da die Kopplungen denen der Neutrinos sehr ähnlich sind und Neutrinobremsstrahlung 
in |HPR,01j schon gerechnet wurde, die nukleonische Seite also schon bekannt gewesen ist, 
ist es notwendig, durch eine Fierztransformation Quark- und Neutralinostrom voneinander 



zu trennen. Friman und Maxwell haben in |FM79|j für die Streuung einen Pion- Austausch 



zugrunde gelegt. C. Hanhart et al. beziehen dagegen die Produktionsrate weicher Strahlung 
(die Neutralinostrahlung wird ja näherungsweise als weich behandelt) auf die on-shell 
A^A^- Streuamplitude (s. |HPR01j ). Dazu wird die Nukleon-Nukleon-Streu-T-Matrix T^]^ 
in Partialwellen zerlegt. Die Partialwellen können in Beziehung zu den Phasenverschiebungen 
aus den Nukleon-Nukleon-Streudaten gesetzt werden. In dem betrachteten Fall ist nur 
elastische Streuung kinematisch erlaubt. 

Für die Berechnung des Neutralinotensors können die Ergebnisse aus dem letzten 
Abschnitt teilweise übernommen werden. 

Durch die Fierztransformation zerfallen die Betragsquadrate der Matrixelemente in 
einen hadronischen Tensor und einen Neutralinotensor: M. = Ti^^M^y. Der Vektorstrom 
trägt in führender Ordnung nicht bei, weil 7° — > 1 und 7 — > in nichtrelativistischer 
Näherung übergehen und der Kommutator mit dem Identitätsoperator verschwindet. 
Denn die Neutralinos können vor oder nach der Streuung der Neutronen abgestrahlt 
werden. In dem einen Fall ist die Energie des Zwischenzustandes positiv, im anderen Fall 
negativ. Die Summe beider Möglichkeiten ergibt den Kommutator zwischen der Streu-T- 




Matrix und dem Operator der Kopplungsstruktur. Der einzige Beitrag stammt von den 
räumlichen Komponenten des axialen Stroms. Demzufolge sind auch nur die räumlichen 
Komponenten des Neutralinotensors von Interesse. Da der hadronische Tensor von C. 
Hanhart et al. in |HPR01j schon berechnet worden ist, beschränke ich mich auf den 
Neut r alinot ensor . 



Von den Quarks zu Neutronen 



Um aus der Kopplung der Quarks an die Neutralinos die Kopplung der Neutronen an 
letztere zu gewinnen, benutzen wir als einfachstes Modell die S't/(6)-Wellenfunktionen 
(s. |Nac86j ) und als Normierung den Neutralen Strom, so daß Renormierungseffekte der 
Starken Wechselwirkung berücksichtigt werden. Dies ist nur eine grobe Näherung und 
vernachlässigt die Dynamik der Quarks innerhalb des Nukleons. Da es sich hier um eine 
Größenordnungsabschätzung handelt, ist dieses Modell aber gerechtfertigt. 



V^qqvVjNC 



46 



KAPITEL 3. NEUTRALINOKÜHLUNG DER SUPERNOVA 



Ts Y 




ur I 
dR -| 



Tabelle 3.3: Isospin und Hyperiadung von up- und down-Quark 

Der effektive Vertex der ggxx-Kopplung sieht nun so aus: 

W^.5xx = ^jMll®YPL-rl(^YPR) (3.35) 

q = k{YI + A' + 2t,YlA); tI^k'Y^ (3.36) 

^ = AT' - N' tan övi. A = cot 6»^/ ^^^-^^ . , ,f 3.37 

T : 3-Komponente des schwachen Isospins (3.38) 

Der Erwartungswert der Hyperiadung wird durch den schwachen Isospin und die Identität 
ausgedrückt: 

If = 57-3; If = 0; 1-1 = ^1, i'H^U'-a + j) (3-39) 

^« - K^ + 5 + H4(5 + -) P.40) 

- '^%^\^) (3-42) 
Die 5'C/(6)-Wellenfunktion des Neutrons ist gegeben durch 

n(l/2,l/2) oc -^i(l{SSv> ^Su^S ^u^SS) (3.43) 
vis 

-[Sd^u' + d^d^v} + d^M^rf^ + d^v^S + w^d^d^ + w^d^ci^)) (3.44) 
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In nichtrelativistischer Näherung geht 7^75 in Gl. ()3.35p aus dem Quarkstrom in er über 
(s. El- Gemäß dem Wigner-Eckart-Theorem braucht man nur den Fall (T3 auszuwerten: 



{n\a3\n) = 1, {nlr^asln) = -\ 



(3.45) 



und dies in (|3.3ö|) einsetzen: 



8M| 



8M| 



K 



K 



1 

3"6 
1 

3"6 



9 
5 
9 



— + - ) 7'^Pi 



(3.46) 



Die Normierung berechnet man analog zu 



24 cos2 Ow sin^ OwMlo 



(3.47) 



weil der Vertex die Struktur r^a^ ® 7^-Pl hat. Einsetzen der letzten Resultate in die 
Ausgangsgleichung ()3.34|) führt zu 



W, 



CaG 



A^F 



nnxx 



i«:^in2övKCos^ÖH.(^y (^ + 3A2 + ^A)yPz. • (3.48) 



Zur besseren Übersicht definiert man einen Faktor ß, der die Kopplungskonstanten enthält: 



V2 



K sin 6w cos 9 



5 '^y-Äü) \T2^'^'^\^ 

Bildung des Betragsquadrats und Kontraktion mit den Spinoren liefert 

4/1 r \ 2 



(3.49) 



8 25 



nHmHwcos^öw{^^ (^^ + 3A2-^A^ (^lü^s + M^)^^^- (3.50) 



wobei berücksichtigt worden ist, daß die linear auftretenden Impulse bei der Winkelintegration 
verschwinden, so daß nur die Energie- und die Massenterme beitragen. Die Z°-Abstrahlung 
wird hier vernachlässigt. 
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3.2.3 Die Emissivität 



Der Energieverlust in ein Phasenraumelement ist gegeben durch 



uJi, UJ2 '■ Energie der abgestrahlten Teilchen 
kl, kl : 3-Impulse der abgestrahlten Teilchen . 




Wenn man die beiden Diagramme aus Tab . 13 . 41 auf summiert . so erhält man ^{[T^n, Ö]), 
da in dem einem Fall die Energie vorher, im anderen Fall nachher abgestrahlt wird. Daraus 
folgt, daß der Identitätsoperator, der ja aus dem Vektorstromanteil der Kopplung stammt, 
nicht beiträgt. 
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In Gl. fl3.51|) fügen wir 1 = / duj6{uj — ui — UJ2) und erhalten: 

-dci;idc<j25(ciJ — cui — üJ2)kik2n{üJi,uj2)^H^-' (3.53) 



nK,.;^) = ^^/3(ü;iü;2 + M3 (3.54) 
Hi = |(M,,p'|[5„T;viv]|M,p)|' (3.55) 

1 2 

Austauschdiagramme sind mit einem Faktor 4 berücksichtigt. Die totale Emissivität eines 
2 — ^ 4- Prozesses hat die Form 

6 = /d^/ ( n ) ^/(^i)/(^2)(l - /(i?0)(l - /(^D) X (3.57) 



X (27r)^5(^) (pi + p2 - - - p'2) ^ (3.58) 

f{E) = p3 , S = - statistischer Faktor (3.59) 

exp(^T^) + 1 8 

Weiterhin wird der 3er-Impuls der Strahlung k gegenüber dem Impuls der Nukleonen in 
den 5-Distributionen vernachlässigt: 

. ^2 

Pn oc y2M]\[T, aus der nichtrelativistischen Energie E = = T (3.60) 

E ^ (X T; die Neutralinos sind relativistisch (3.61) 

^ = ^ 0.12 mit T^30MeV, MTv^ilOOOMeV . (3.62) 

PN V2M^ 

Dies vereinfacht die Integrationen, die dann nur noch den Neutralinotensor betreffen. 
Dabei wird die sphärische Symmetrie und Energie-Impuls- Erhaltung ausgenutzt. 

Im nächsten Schritt definieren wir Gesamt- und Relativimpulse sowie dimensionslose 
Variablen: 

1 1 
P = Pi + P2 = pi + P2, P = ^{Pi - P2), P' = -(pi - P2), (3.63) 

2 P' /2 P" ,2 /i-M y 

2MT' 2MT' 8MT' ^ t ' T ^ ' 
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T ist die Temperatur des Gores und damit die Temperatur des Neutronengases, M 
die Masse eines Neutrons. Die z-Achse wird entlang von P gewählt. Außerdem können 
einige Winkelintcgrationcn schon ausgeführt werden. Bei hohen Temperaturen wird die 
Winkclabhängigkcit in den Fcrmi-Funktionen schwach, im Grcnzfall völliger Entartung 
verschwindet diese Abhängigkeit. Da das Neutronengas nahezu entartet ist, können die 
Fermi-Punktionen durch ihre winkelgemittelten Funktionen ersetzt werden: 



dcos7 



-1 



+ 1 



^u^+A^— 2uAcos7 



+ 1 



In 



2Xusmh{u^ + X^-y) \^cosh - A)2 - y)) 
Genauso wird ein winkelgemittelter hadronischer Tensor eingeführt: 



(3.65) 
(3.66) 



7ii 



(3.67) 



diese Integrationen können analytisch ausgeführt werden, da die T-Matrix als Partialwellenzerlegun^ 
vorliegt. Übrig bleibt eine vierdimensionale Integration: 

CO oo 

e = ß^^{2M^Y-'T'-'M'^ j d6e-'h{5/M^) J daF{6, a)n{Ta)u, (3.68) 

mit den Abkürzungen 

6 

h{x) 
h{x) 

m <t) 

L(a, X) 



p2 _ pl2 p2 pl2 

2MT ' ^ 



X 



2MT ' 

[{x''-l)Ioix)-ix-lfl2{x)], 



-1 

JdXL y^{a + S),X^ L y^{a-S),X 



cosh(|((A + a)2 -y)) 
sinh(a2 + X^-y) Vcosh(|((A + a)2 - y)) 



In 



(3.69) 
(3.70) 

(3.71) 

(3.72) 
(3.73) 
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10° 



10^ 



1 10 



T3 



5 10" 



2-10-^ 




10" I — ^ — ^ — ^ — ^ — ^ ^ ^ ' — ^ 

2 4 6 8 10 

Abbildung 3.8: Emissivität der Bremsstrahlung 

Die Funktionen /^(x) lassen sich durch Taylorentwicklung der Wurzelterme bis zur dritten 
Ordnung in {t/xY recht genau analytisch berechnen. Die Entwicklung hat sogar bei x = 1 
nur einen Fehler von 5% auf das exakte Ergebnis. Alle anderen Integrationen werden 
numerisch mittels verallgemeinerter Gauß- Integration ausgeführt. Als Entartung r] folgt 
für T = 30 MeV und Kerndichte ein Wert von 1.3. 

Das Endergebnis 

Das Endergebnis läßt sich mit hoher Genauigkeit durch 

e = a(1.5 + 6M^ + 8.6M^)e-^^'^>^ x 10^^-^, Tcore = 30 MeV (3.74) 

^ g-s 

approximieren (s. Abb. üTHI) . 

Wegen der soft-radiation-Hah.eixm.g ist dieses Ergebnis für alle Felder gültig, die an 
den axialen Strom koppeln. Man muß nur das entsprechende ß bestimmen. 

Massenschranke für ein Bino 

Mit diesem Resultat kann eine untere Massenschranke für ein Bino in Abhängigkeit der 
Squarkmasse berechnet werden, siehe Tab. 13.51 

Dieses Ergebnis widerspricht nicht der Existenz eines leichten Neutralinos mit einer 
Masse von 34 MeV, wie es z. B. in |(]DR,SOÖ] vorgeschlagen wird, wenn die Squarkmasse in 
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Temperatur 


Squarkmasse 


Massenschranke 


(MeV) 


(GeV) 


(MeV) 


30 


100 


147 


30 


200 


110 


30 


500 
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Tabelle 3.5: Untere Massenschranke 



der Größenordnung von einigen TeV liegt. Im Vergleich zu den unteren Massenschranken, 
die man aus dem Elektron- Positron-Kanal erhält, sind die letzteren etwas kleiner, liegen 
aber in derselben Größenordnung. 



Kapitel 4 
Diskussion 



In diesem Kapitel möchte ich die Ergebnisse und die in sie einfließenden Annahmen aus 
dem vorhergehenden Abschnitt eingehender untersuchen. Die beiden Verfahren sollen 
miteinander verglichen werden, und die Neutralinoemissivität soll der Neutrinoemissivität 
gegenüber gestellt werden. 



4.1 Diskussion des Elektron-Positron-Kanals 

In die Berechnung der Emissivität des Elektron-Positron-Kanals sind eine Reihe von 
Eingangsdaten eingegangen, die im folgenden noch etwas ausführlicher diskutiert werden: 

1. Temperatur- und Entartungsprofil einer Supernova 

2. Coredichte 



3. Dauer der Neutralinoemission und Menge der abgestrahlten Energie 

4. Selektronmasse 

5. Mischungsparameter 

6. Näherungen 



Temperatur- und Entartungsprofil einer Supernova 

Das Temperatur- und Entartungsprofil einer Supernova stammt aus |BL86j . Die Autoren 
haben eine Supernova simuliert, dazu mußten sie Anfangs- und Randwerte annehmen, um 
die Differentialgleichungen zu lösen. Uber die Anfangswerte einer Supernova-Explosion 
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weiß man wenig. Kleine Änderungen der Anfangsdaten führen zur gleichen Entwicklung 

nach etwa einer halben Sekunde. Die Supernova befindet sich nach dieser Zeit im Diffusionsgleichgewic 

Als Temperatur- und Entartungsprofil habe ich die Profile nach einer halben Sekunde 

genommen. Man kann somit vermuten, daß dies vernünftige Profile sind. Für meine 

Rechnungen habe ich die Profile über die Zeit als konstant angenommen (es liegen auch 

nicht mehr Daten im relevanten Zeitintervall vor), betrachtet man deren Abbildungen, so 

stellt man fest, daß die Änderungen nach einer weiteren halben Sekunde nicht allzu groß 

sind, vor einer halben Sekunde allerdings ist die Änderung, aber auch die Unsicherheit 

groß. Diese Profile sind ungefähre Mittelwerte für die erste Sekunde. 

Coredichte 

Da die Verteilungen in Abhängigkeit der eingeschlossenen Baryonenmasse angegeben sind, 
muß die Volumenintegration in eine Massenintegration transformiert werden. Dafür wird 
die Coredichte benötigt. Da das Gore Kerndichte erreicht hat und diese Dichte bekannt 
ist, stellt dies kein Problem dar. Die Dichte wird konstant gehalten, was nur nicht ganz 
richtig ist. Weil aber die Emissivität aus einer dünnen Kugelschale stammt, macht das 
nichts aus. 

Dauer der Neutralinoemission und Menge der abgestrahlten Energie 

Die abgestrahlte Gesamtenergie kann einfach aus einer kollabierenden, kugelsymmetrischen 
Massenverteilung abgeschätzt werden. Das meiste davon wird in Form von Neutrinos 
abgestrahlt. Ein Teil allerdings soll in Form von Neutralinos abgeführt werden. Da anzunehmen 
ist, daß die Neutralinos nicht wie die Neutrinos gefangen werden, können erstere die 
Supernova sofort verlassen. Wir wissen aber, daß der Hauptteil der Energie von Neutrinos 
abgeführt wird. Also dürfen die Neutralinos nur eine kleine, im Verhältnis zu den Neutrinos 
vernachlässigbare Störung darstellen. Die Neutrinos führen einige 10^^ erg ab, die Neutralinos 
sollen maximal etwa 10^^ erg abführen. Da die Unsicherheit bezüglich dieser Zahl groß ist, 
habe ich verschiedene Zeiten gewählt, um aus der Leistung P = J dMsteme die Energie 
E = PAt zu berechnen. Vergleicht man dies mit dem Kriterium von Raffelt ( 'Raf 96J ) , der 
eine Höchstemissivität von 10^^erg/(gs) für zulässig hält, ohne daß das Neutrinosignal 
gestört wird, so ist meine Schranke um einen Faktor Drei niedriger (10^^ erg/(gs) ■ 2, 8 • 
XO^^g .13 = 2.8 ■ 10^^ erg). Bei diesem Vergleich muß man vorsichtig sein, weil dieses 
Kriterium auf diesen Prozeß nicht anwendbar ist. Die Neutralinos werden im Gegensatz zu 
den Neutrinos nicht verzögert emittiert. Dies ist aber eine Voraussetzung zur Anwendbarkeit 
des Raffelt-Kriteriums. Da ich von einer Emissionsdauer von etwa einer Sekunde ausgehe, 
habe ich die Profile nach einer halben Sekunde gewählt, so daß diese Profile ein zeitliches 
Mittel darstellen. 
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Selektronmasse 

Uber die Selektronmasse wissen wir nur, daß sie mindestens im 100 GeV-Bereich liegen 
muß. Ich habe drei verschiedene Werte für die Selektronmasse gewählt und eine Massenschranke 
berechnet. Wird die Selektronenmasse zu groß, wird die Berechnung einer unteren Massenschranke 
nicht mehr möglich. Ein Vorteil der Rechnung ist, daß die Selektronmasse und die Mischungsparameter 
(s. u.) in einem Faktor stecken, der nicht in die Integration eingeht. 

Wenn die Selektronmasse im TeV-Bereich liegt, steht ein leichtes Neutralino (30MeV, 
ICDRSOO] ) nicht im Widerspruch zu Supernova- Ereignissen. 

Mischungsparameter 

Ebenso unbekannt sind die Mischungsparameter Nij (Mischungsmatrix) der Neutralinos. 
Aus den Experimenten ergeben sich aber Schranken an diese Parameter (siehe vorheriges 
Kapitel). Aus der Supernova lassen sich keine Folgerungen darüber ziehen. Der Einfluß 
auf die Kopplung ist gering, solange die Sfermionmassen im Bereich der Z°-Masse liegen. 
Bei Sfermionmassen im TeV-Bereich wird der Einfluß der Ni^ß größer. 

Näherungen 

Bei der Berechnung des Wirkungsquerschnitts sind eine Reihe von Näherungen gemacht 
worden: Die Elektronmasse wurde zu Null gesetzt, die Impulsüberträge in den Propagatoren 
vernachlässigt und die Massen von Selektron(L) und Selektron(R) als gleich angenommen. 
Die ersten beiden Näherungen sind durch die Kleinheit der Größen gerechtfertigt. Die 
letzte Näherung ist eine Vereinfachung, die mit der Supergravitation gerechtfertigt werden 
kann. 

4.2 Die Nukleon-Nukleon-Bremsstrahlung 

Einige Annahmen, die für den Elektron-Positron-Kanal gemacht werden müssen, fallen 
hier weg, da bei der Nukleon-Nukleon-Bremsstrahlung Emissivitäten mit dem Raffelt- 
Kriterium verglichen werden. Die Partialwellen der Nukleon-Nukleon-Streuung sind dem 
Experiment entnommen. Die Streuung wird realistischer behandelt. Für die Näherungen 
gilt dasselbe wie für den Elektron-Positron-Kanal. 

Die Unbekannten stecken in dem Faktor ß: die Mischungsparameter und die Squarkmasse 
(der Z*^-Kanal wurde in dieser Rechnung vernachlässigt), ß verhält sich wie 




(4.1) 
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(Die Z^-Massenabhängigkeit stammt aus der Normierung durch den Neutralen Strom). 
Auch hier ist die wesenthche Abhängigkeit durch die Squarkmasse gegeben. Ein leichtes 
(auch masseloses) Neutralino ist möglich, wenn die Squarkmasse im TeV-Bereich liegt. 

Soft-radiation-Näherung 

Die abgestrahlten Neutralinos wurden als soft behandelt. Dies ist schon für leichte Neutralinos 
(cä 30 MeV) problematisch. Die Nukleon-Nukleon-Streuung spielt sich auf einer Längenskala 
ab, die einer Pionmasse (~ 100 MeV) entspricht. Der Entwicklungsparameter der Näherung 
ist also mit ~ 4 schon für ein leichtes Neutralino nicht besonders klein. Damit ist diese 
Näherung nicht gerechtfertigt. Aber es gibt nichts Besseres. Außerdem wurden Effekte 
vernachlässigt, die durch das Medium Kernmaterie verursacht werden. 

4.3 Vergleich mit der Neutrinoluminosität 

In Abb. 2 aus jBurQOj ist die Neutrinoluminosität in den ersten Millisekunden nach 
dem Kollaps dargestellt. Der Peak bei 113 ms stammt von den durch die Schockwelle 
ausgetriebenen Elektron-Neutrinos. Die Luminosität erreicht 6 ■ lO^^erg/s, die im Peak 
abgestrahlte Energie beträgt aber nur 2-10^^ erg, das Spektrum ist ein Einfangspektrum, 
kein (pseudo-)thermisches und damit auch kein Schwarzer Strahler. Zu derselben Zeit setzt 
abrupt die Emission von Myon-Neutrinos und Elektron-Antineutrinos ein, die kleinere 
Luminositäten erreichen. Nach dem Durchlaufen der Schockwelle ist das Spektrum ein 
pseudothermisches. 

Da das Gore für Neutralinos durchsichtig ist, kann man solche Peaks in der Luminosität 
für die Neutralinos nicht erwarten. Die Emissivität wird vielmehr durch das Zusammenspiel 
von Gore- Temperatur und dem chemischen Potential von Elektronen und Positronen 
getrieben. Die Emissivität findet in meinem Modell einer Schale bei einer eingeschlossenen 
Masse von einer Sonnenmasse statt. An der Stelle ist die Temperatur hoch und das 
chemische Potential der Elektronen niedrig. 

In |KacOOj wurde die Neutralinoemissivität bezogen auf das Gorezentrum für drei 
verschiedene Entartungen berechnet (das hat den Vorteil, daß Mischungsparameter, Temperatur, 
Selektronmasse nicht gebraucht werden). Für ein Neutralino mit einer Masse von 34 MeV 
ist die Emissivität kaum unterdrückt. Kachelrieß fordert in jKacOOj eine deutliche Unterdrückung 
(M^/T > 10), folglich muß entweder die Temperatur niedrig sein oder die Masse des 
Neutralino hoch. 

Anmerkung Abb. 1 in [KacOOj konnte von mir nicht exakt reproduziert werden. Für 
größere x = M^/T- Werte habe ich höhere relative Emissivitäten erhalten: Bei x = 2 



4.3. VERGLEICH MIT DER NEUTRINOLUMINOSITÄT 
ergibt sich R — e{x)/e{0) ^ 0.7 statt R 0.6 wie bei Kachelrieß. 



Anhang A 



Der totale Wirkungsquerschnitt 

a(e+e" XiXi) 

In diesem Abschnitt soll die Berechnung des totalen Wirkungschnitt s des Prozesses e~^e~ 
XiXi dargestellt werden. Dabei stütze ich mich auf die allgemeinere Rechnung von A. Bartl 
et al. (IHFMMI). 

A.l Beteiligte Diagramme 

Folgende Prozesse tragen in niedrigster Ordnung zum Wirkungsquerschnitt bei: 




A.2 Feynmamegeln 

Die Feynmanregeln für die Vertizes des Standardmodells sind "HM84I entnommen, bei den 
Regeln für die SUSY- Teilchen beziehe ich mich auf ^HK85j und |.GH86j . Die Elektronmassen 
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werden vernachlässigt sowie die Impulsüberträge gegenüber den Massen des Z'^-Bosons 
und des Selektrons. Weiter sollen links-chirales und rechts-chirales Selektron (d. h. das 
dazugehörige Standardmodellelektron ist links- bzw. rechts-chiral) in ihren Massen entartet 
sein. Als Neutralinobasis wird die Photinobasis verwendet. 



Die Propagatoren sehen wie folgt aus {i: Fermion mit schwachem Isospin —1/2): 

i i 



-Propagator 



Z -Propagator 



M| 



p2 - M| 



p2 - M2 



ZO 



(A.l) 
(A.2) 



Für den Z^-Propagator wurde die Landau-Eichung benutzt. Die Vertizes stehen in Tabelle 
lA.ll Die dabei auftretenden Konstanten bedeuten: 

^(cy + Ca); Cr = ^{cv- Ca); cy = T3- 2e£sin2 6w; ca = T3 (A.3) 



cl 
Li 



\ [N^.Nl, - N,,N*,,) cos(2/3) - i {N,,N*, - N,,N*,) sin(2/3) = -R, . (A.4) 



Zur Abkürzung vereinbaren wir: 



h 



ri 



V2 
V2 



] sin öwNii + 



1 



cos 9w 
; sin 9wN^i — Ci 



w 



1 
2 

Nt 



Ce sin^ 9w 



N- 



12 



cos 9w 



(A.5) 
(A.6) 



• • 

A.3 Die Ubergangsmatrixelemente 

Die Anwendung der Feynmanregeln liefert fünf Ubergangsmatrixelemente: 

-iMl = -^^Dzo{s)v{pi)j''{cLPL + CRPR)u{p2)ü{p3)jf,{LiPL + RiPr)v{p4), 



(A.7) 

-iM2 = -9^\li\^D,{t)v{pi)PRv{p3)ü{p4)PLu{p2), (A.8) 
-iMs = -g%\^D,{u)v{pi)PRv{p^)ü{p3)PLu{P2) ■{-!), (A.9) 
-iMi = -gyi?D,{t)v{pi)PLv{p3)u{pA)PRu{p2), (A.IO) 
-iM^ = -gy^\^D,{u)v{p^)PLv{p^)u{p3)PRu{p2) ■{-!), (A.ll) 

mit Dz.{s) = -^, D,{x) = -^ . (A.12) 



ANHANG A. DER TOTALE WIRKUNGSQUERSCHNITT a{e+e- x 



Graph des Vertex 



Vertexfaktor 



X? 




X? 




X? 




xS 




e cos 



Xi Xi 



cos 



Tabelle A.l: Vertexfaktoren 
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Da die Ausgangszustände identische Fermionen sind (Neutralinos sind Majorana- Fermionen), 
erhalten die Amphtuden für M3 und M5 das zusätzhche, angeschriebene Minuszeichen. 



A.4 Der differentielle Wirkungsquerschnitt 



A.4.1 Berechnung der Betragsquadrate 



Der nächste Schritt ist die Berechnung von 




(A.13) 



(T(out) 



CT (out) 




(A.14) 
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Hierbei bedeutet (a-b) = a^J)^. Die Summation über die auslaufenden Spins wird ab sofort 
nicht mehr mit geführt und imphzit vorausgesetzt. 

•^1-^^ = 7Tr^£rZ- ciL2(pi.p4)(p2-p3)+4^l(pi-P3)(p2-p4), (A.15) 
M^o COS \ 

+ clL^RiM^p^ ■ P2) + clL^RiM\p^ ■ p^) 

+ cMd'l ■1'4)(P2 -Ps) + cIRKpi ■P3)ip2-P4)) 

A.16) 

A.17) 
A.18) 
A.19) 

A.20) 

A.21) 

A.22) 
A.23) 

A.24) 

A.25) 
A.26) 

A.27) 
A.28) 
A.29) 



cos^ 9w 

1| . |2 , 

MlMt = - ,.2 V 2ü ('^crRiM\p,-P2)+8crLi{p,-P3){p2-P4)), 
IVI ^QlVIg2 cos n]v \ / 



■>4|; 14 



■M2Atr = ^J^4(pi-p3)(P2-P4), 



e 



A^2A^| = 0, 
M2Mt = 0, 



A^3A^+ = 0, 
MaMt = 0, 

M^Mt = ^^4(pi-p3)(P2-P4), 

e 
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p3 



pi 





p2 



p4 

Abbildung A.l: Impulse im Schwerpunktsystem 

A.4.2 Weitere Vollständigkeitsrelationen 

Für die Berechnung des Betragsquadrates des Gesamtmatrixelements sind die folgenden 
Vollständigkeitsrelationen nützlich. 

"(P' ^)^iP^ s)=f + m, ^ t;(p, s)tJ(p, s)=f-m, (A.30) 

s s 

^^;^(p,s)fZ(p,s) = C-i(|^ + m), J]T;(p,s)ii^(p,s) = (|^-m)C^ . (A.32) 

s s 

Diese Regeln werden sofort einsichtig, wenn man 

u{p,s) = Cv^{p,s); v{p,s) = C'iF{p,s) (A.33) 
in die bekannten Vollständigkeitsrelationen einsetzt. Eine weitere hilfreiche Identität ist 

= -C'-^'^C . (A.34) 

A.4.3 Die Kinematik 

Der totale Wirkungsquerschnitt ist ein Lorentzskalar. Also kann man ein beliebiges Bezugssystem 
wählen, um den totalen Wirkungsquerschnitt zu berechnen. Ich wähle als Bezugssystem 
das Schwerpunktsystem, weil dann die Rechnung einfacher ist. 
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XiXi 



Als Parametrisierung der Impulse (genauer: 4er-Impulse) wähle ich (siehe Abbildung 

EU 

Pi={E, 0, 0, E,) ,P3={uJ, 0, ksinO, kcosO,) , (A.35) 
P2 = {E, 0, 0, -E,) ,P4=(^, 0, -ksinO, -kcosO,) . (A.36) 

Hierbei sind k = \p^\ = \p^\ der Betrag des 3er-Impulses, uj die Energie der auslaufenden 
Teilchen, E die Energie und der Impuls des einlaufenden Elektrons und Positrons (Sie 
werden als masselos betrachtet.). Damit folgt für die Mandelstam- Variablen: 

s = {pi+P2? = ^E^ , (A.37) 

t = (pi - p3)2 = _i (^s _ ^s{s - 4M2) cos + , (A.38) 

(j9i - = -1 (^s + ^s{s - 4M2) cos + . 



u 



(A.39) 



Die Minkowski-Produkte in den Gleichungen ()A.16|) - ()A.29|) werden mit diesen Beziehungen 
durch die Mandelstam- Variablen ausgedrückt: 



{Pl 


■P2) 


1 

= 2^ ' 




(A.40) 


{Pl 


Ps) 


= {.P2-Pa) = -i^{t- 




(A.41) 


{Pl 


■Pi) 


= {P2-P3) = -^{U 




(A.42) 


{P3 


■Pi) 


= -s-M^ . 
2 




(A.43) 



A.4.4 Der difFerentielle Wirkungsquerschnitt 



Der differentielle Wirkungsquerschnitt im Schwerpunktsystem für eine 2 
ist allgemein gegeben durch (siehe z. B. [IIM84J): 

der 1 Pf 
dQ GAtc'^s Pi 



2-Streuung 
(A.44) 



Außerdem muß noch über die Eingangsspins gemittelt werden. Da die auslaufenden Neutralinos 
identisch sind, ist der Phasenraum nur halb so groß. Einsetzen liefert: 



da 

dn 



9^ 



1287r2s 

2 t2 



4M2 



(t - M^)2 + {u- M^f - 2M^s 



X 



4^1 + 4^1 Li{(^l\i 

M|o cos4 e 



w 



Ml, Ml cos2 Ow 



AMi 



(A.45) 
(A.46) 
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A.4.5 Der totale Wirkungsquerschnitt 

Der totale Wirkungsquerschnitt ist gegeben durch o" = /q ^ dfi. Weil wir es in unserem 
Falle mit 2 — *■ 2-Streuung zu tun haben, hängt der differentielle Wirkungsquerschnitt nur 
vom Winkel zwischen einlaufendem und auslaufendem Teilchen ab. Für die Integration 



ist nur der Term {t-M^y + 



[u- 



■M' 



2\2 



2M^s relevant, die restlichen Terme sind für die 



Integration unwichtige Vorfaktoren. Setzt man die Gleichungen ()A.38jl und ()A.39|1 für die 
Mandelst am- Variablen ein, so vereinfacht sich der differentielle Wirkungsquerschnitt zu 



da 



cos 



1287r2 



M|o cos4 Ow MloM? cos2 öw 



\h\' + \ri\' 
AMI 



X 



3/2 



(A.47) 



Der Wirkungsquerschnitt a hat dann folgende Gestalt, wenn man /^(H-cos^ 6')d0 dcos 9 



yvr benutzt: 



3 sin"^ 6w 



Ml, Ml cos2 Ow 



'2 

M|o cos4 Ow 



+ 



1 - 



4M| 
4m^ 



3/2 



S. 



(A.48) 



Anhang B 

Fierz-Transformation 



B.l Fierz-Transformation 

B.1.1 Allgemeine Transformationsformel 

Ziel dieser Transformation ist es, Strom- Strom- Wechselwirkungen der Form 

Mab oc {viTaV2) {Ü3TBU4), Tk e {1, 75, 7m> 757m> ^-^i- = ^/Hlßli^ - 71^7^)} 
wie sie zum Beispiel in der Reaktion e'^e~ Vv auftreten, in der Form 



zu schreiben (Die Summe läuft über alle Lorentz-Kovarianten). Dann zerfallen auch 
die Interferenzterme mehrerer Matrixelemente in Produkte zweier Tensoren mit gleicher 
Stromstruktur, was sonst nicht immer der Fall ist. Dadurch können Rechnungen einfacher 
werden. Man entwickelt {vi Ta V2) (ös «4) nach einem vollständigen Satz von Matrizen, 
nämlich den 16 F-Matrizen. Außerdem benutzt man, daß diese orthogonal zueinander 
bzgl. der Spur als Skalarprodukt und (anti-)hermitesch sind. 
Die Entwicklungskoeffizienten ergeben sich wie folgt: 




(B.2) 



C,D 




(B.3) 
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Fierz- Transformation 




Tabelle B.l: Wirkung der Fierz- Transformation 



Die Gramsche Matrix Gab — {T^(^a^b)) lautet (Reihenfolge: Skalar, Pseudoskalar, 
Vektor, Axialvektor, Tensor): 



(l 
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Vo 
















(B.4) 



B.l. 2 Anwendung auf qq — ^ 



Die Fierz- Transformation wird auf das Matrixelement M. = u{pi)PRv{j)4)u{p^)Pi,u{p2) 
angewandt, was der Struktur |(1 + 75)®|(1— 75) entspricht. Es ist 



^(1+75)® ^(1-75) 



-(1 ® 1 - 1 ® 75 + 75 ® 1 - 75 ® 75) 



(B.5) 
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ANHANG B. FIERZ-TRANSFORMATION 



Anwendung von Formel IB. 31 auf die einzelnen Terme ergibt: 

1 ® 1 = ^ (^1 ® 1 + 75 ® 75 + 7m ® 7^ - 757m ® 757^" + ^c^ßu ® cr^""^ , (B.6) 

1®75 = ^ (^1®75 + 7m®757'' + 75® 1 -757m®7'' ^ , (B.7) 

75 ® 1 = ^ (^1 ® 75 - 7m ® 757^* + 75 ® 1 + 757m ^ 1^ ^ , (B.8) 

75 ® 75 = ^ (^1 ® 1 + 75 ® 75 - 7m ® 7'' + 757m ® 757^" + ^c^ßu ® (r^""^ ■ (B.9) 
Man erhält nach Einsetzen dieses Zwischenresultats in Gl. (IHI 



® = ^(1 + 75) ® ^(1 - 75) = ^7m(1 + 75) ® 7^1 - 75) = \i^.Pr ® YPl (B.IO) 

Dieses Resultat ist aufgrund der Symmetrieeigenschaften der Ausgangsstruktur nicht 
weiter verwunderlich! 

B.2 Nichtrelativistischer Grenzfall der Kopplungsstrukturei 

Ein Diracspinor M*^*^(p), > 0, oder < 0, hat als Lösung der Dirac-Gleichung 

in der Standard-Darstellung die allgemeine Form 

= (J) ' = (?) • ^^-12) 

Wenn die kinetische Energie gegenüber der Masse vernachlässigbar klein ist, sind die 
kleinen Komponenten der Diracspinoren unterdrückt: 



E+m 



Dadurch überleben von den fünf Lorentzkovarianten nur solche, die große mit großen 
Komponenten verbinden. Am Ende bleiben nur zwei Kopplungsstrukturen übrig. Dazu 
verwendet man die 7-Matrizen in der Dirac-Darstellung: 



_ M 0\ - ^\ 5 / ^ 



B.2. NICHTRELATIVISTISCHER GRENZFALL DER K0PPLUNGSSTRUKTURENQ9 



^fc^5^ ^ij gjj^j-^ diagonal, alle übrigen Kombinationen mischen kleine mit großen Komponenten 
und verschwinden deshalb (vgl. dazu ICXilQSj ). 



üu - 


- x' 


'x, 






(B.15) 


üj^u - 




^x, ^ 




0; sonst, 


(B.16) 


üj^u - 


- 0, 








(B.17) 


ü'j^'y^u - 


- x' 






= i] sonst. 


(B.18) 


üa^^u - 


- x' 






= j, u = k,ijk, zyklisch; sonst. 


(B.19) 



(Griechische Buchstaben bezeichnen Indizes, die von bis 3 laufen, römische Buchstaben 
Indizes von 1 bis 3, m ist ein Diracspinor, x ein Weylspinor.) Man sieht, daß nur / und er 
übrigbleiben. 



Anhang C 

Partialwellenzerlegung 



Hier wird die Partialwellenzerlegung des winkelgemittelten hadronischen Tensors Ti 
angegeben. 



{SM'sp'\Tr,N\SMsp) 
T{p,S'M's,SMs,L',AL,Mj) 

S,Ms 
L, Ml 
J,Mj 
hvLS 



Yl yL'M',{p)YLM,{prnp, s'm's, sms, al, mj) 

(C.l) 

^{S'M's, L'M'l\JMj){S'M's, L'M'l\JMj)T^^'^%p) 

(C.2) 

^^""^^ (^) ^'''""^'^ {hvLsip)) (C.3) 
(falls die Kanäle ungekoppelt sind) 
Spin und dazugehörende 3-Komponente 
Drehimpuls und dazugehörende 3-Komponente 
Gesamtdrehimpuls und dazugehörende 3-Komponente 
Phasenverschiebung der Partialwelle JL'LS . 



Eine einfache, aber längliche Rechnung führt zur Beziehung zwischen Ti, und T: 



(C.4) 



-29^e(T(p, l{M's - 1), IM5, L', AL, Mj - l)*T(p, IM^, l(Ms + 1), L', AL', Mj) ) 
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Anhang D 

Liste der verwendeten Symbole 
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ANHANG D. LISTE DER VERWENDETEN SYMBOLE 



Symbol 


Bezeichnung 


SI-Einheit 


natürliche Einheit 


s 


Länge 


m 


eV-i 


t 


Zeit 


s 


eV-i 


a 


Bcsciilcunigung 


ms~^ 


eV 


E 


Energie 


J 


eV 


F 


Kraft 


N 


eV 


L 


Leuchtkraft 


kg/s 


eV^ 


M 


Masse 


kg 


eV 


P 


Druck 


Nm-2 


eV^ 


p, k 


Impuls 


Ns 


eV 


R, r 


Radius 


m 


eV-^ 


T 


Temperatur 


K 


eV 


u, e 


Energiedichte 


Jm~^ 


eV^ 


V 


Volumen 


m 


eV-3 


V 


Geschwindigkeit 


ms~^ 


1 


e 


Emissivit ät 


Jm~^s~^ 


eV^ 


P 


Dichte 


kgm~^ 


eV^ 


CT 


Wirkungsquerschnitt 


b 


eV-2 


P 


Leistung 


W 


eV-2 


ß 


chemisches Potential 


J 


eV 


V 


Entartung 


1 


1 


w 


"chemische Zusammensetzung" 


1 


1 


N 


Teilchenzahl 


1 


1 


n 


Teilchendichte 


m~^ 


eV=^ 



Tabelle D.l: Liste physikalischer Größen 
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Syiii1)()1 


Konstante 


AAert 




e 


Elftmftnt ar ladi m et 


1 60 • 10""^^ As 






Ppinstriiktiirkonstantp 


1/137 




G 


Gr avit at ionskonst ant e 


6.67- IQ-" N ■ m^ 


kg-^ 


k 


Boltzmann-Konstante 


1.38- 10^23 j/j^ ^ 


= 0.862 - 10-W/K 


h 


Planck-Konstante 


6.625- 10-3^ Js 




c 


Lichtgeschwindigkeit 


3 • 10« m/s 




Mo 


Sonnenmasse 


1.99 ■ 10^° kg 




m„ 


Masse Neutron 


1.675- 10-27 kg = 


939.6 MeV 


rup 


Masse Proton 


1.673- 10-27 kg = 


938.3 MeV 


9w 


schwacher Mischungswinkel 


sin^ Ow = 0.223 





Tabelle D.2: Liste der Konstanten 



Symbol 


Operator 


möghche Darstellung 


C 


Ladungskonjugation 


17270 


Pl,r 


Projektoren 


1(1 T 75) 


a* 


komplexe Konjugation von a 






Transposition von A 






hermitesche Konjugation von A 





Tabelle D.3: Liste der Operatoren 



Symbol 


Bezeichnung 


C 


Lagrangedichte 


M 


Matrixelement 


Mj 


Ncutralinotensor 




Hadrontcnsor 


N 


unitäre Transformation 



Tabelle D.4: sonstige Symbole 



ANHANG D. LISTE DER VERWENDETEN SYMBOLE 



cgs — > natürliche Einheit natürhche — > cgs Einheit 



1 cm 


= 5,068- lOM/eV 


leV = 


5.068 


10^1/cm 


Is 


= 1.52- 10^5 1/eV 
= 8.62- 10-5 eV 


leV = 


1.519 


lO^^l/s 


IK 


leV = 


1.160 


10^ K 


1 erg 


= 6.24- lO^^eV 


lcV = 


1.602 


10"^2 erg 


lg 


= 5.61 • 10^2 eV 


leV = 


1.783 


10-33 g 



Tabelle D.5: Umrechnung cgs-Einheiten-natürliche Einheiten 
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